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ABSTRACT 


This  thesis  demonstrates  a  remote  terminal 
approach  to  computer  mapping,  and  suggests  that  the 
terminal  should  be  utilized  in  future  work  in  computer 
cartography. 

A  package  was  developed  in  APL  primarily  to  help 
landscape  architects  and  urban  planners  in  representing 
combinations  of  various  environmental  factors  by  maps. 


Factors  such  as  rivers,  forests,  and  slopes 
are  represented  in  the  computer  as  matrices  called  base 
maps.  The  X^.th  element  of  each  matrix  corresponds  to 
the  (i,j)th  point  of  a  rectangular  grid.  If  a  factor 
present  at  the  (i,j)th  point,  the  corresponding  elemen 
of  the  area  is  set  to  1  ,  otherwise,  it  is  set  to  0  . 


is 

t 


Various  logical  and  arithmetic  operations  are 
available  to  combine  base  maps  into  new  composite  maps. 

These  show  meaningful  results  that  could  not  be  easily 
produced  by  hand.  All  APL  functions  can  also  be  used. 

They  can  be  applied  simply  and  naturally.  Thus  the 
user  has  a  versatile  and  powerful  tool  for  acting  upon 
his  maps. 

Functions  to  facilitate  base  map  input,  manipulation, 
generation  and  output  were  developed  and  incorporated  in  the 
package.  Functions  to  produce  contour  map  output  are  also 
provided.  All  outputs  use  APL  characters  to  represent  levels 
of  gray.  The  number  of  gray  levels  is  limited  only  by  the 


number  of  different  characters. 
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The  use  of  the  package  can  be  extended  to  other 
similar  problems  and  fields  wherever  the  information  can 
be  displayed  in  the  form  of  a  map. 
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INTRODUCTION 


In  1942,  an  eminent  American  geographer  wrote: 
"Maps  are  drawn  by  men  and  not  turned  out  automatically 
by  machines".  Yet  today,  using  computers  and  computer 
oriented  techniques,  cartography  has  a  new  and  efficient 
tool  for  producing  maps  and  representing  spatially 
distributed  data. 

One  of  the  new  approaches  that  can  be  used  in 
real-life  geography  problems  is  suggested  and  elaborated 
in  this  work. 

The  aim  of  the  work  is  primarily  to  develop  a 
small  package  for  urban  planners,  landscape  architects  and 
designers,  so  that  they  can  enter  and  combine  various 
spatially  distributed  environmental  factors  and  display  the 
result  in  the  form  of  a  map. 

As  the  design  process  involves  interaction  between 
a  designer  and  the  objects  of  his  design,  it  is  intended  not 
to  put  the  computer  between  the  designer  and  these  objects, 
but  to  make  the  computer  an  invisible  communication  link 
and  tool  for  him.  This  will  allow  him  to  act  on  the  objects 
of  his  design  without  breaking  the  design  process  either 
in  space  or  in  time.  He  will  be  able  to  see  the  result 
immediately. 

This  could  be  achieved  by  doing  the  interactive 
technique.  There  are  two  ways  to  do  it,  either  by  using 
a  CRT  or  by  a  remote  typewritten  terminal.  The  2741 
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terminal  was  used  because  it  was  readily  available.  This 
choice  has  proved  satisfactory  for  the  following  reasons: 

1)  the  cost  associated  with  use  and  purchasing  or 
leasing  the  terminal  is  convenient  for  the  user. 

2)  hard  copy  output  suitable  for  xeroxing  is 
obtained  immediately  on  a  remote  terminal. 

3)  similar  work  using  APL  and  remote  terminal  has 
not  been  reported  anywhere  to  the  best  of  my  knowledge. 

4)  the  efficiency  and  suitability  of  APL  is  very 
appealing . 

5)  remote  terminal  can  be  installed  wherever  there 
is  a  telephone  outlet. 

6)  there  is  no  need  for  punch  cards,  paper  tape  or 
similar  intermedia  for  programs  and  data. 

7)  future  communication  networks,  data  banks 
(urban  data  banks),  computer  utilities,  etc.  are  very 
promising  in  view  of  remote  access  via  the  terminal. 

The  intent  of  the  work  is  to  develop  the  package 
to  the  state  where  it  can  be  used  effectively  by  the 
designer  without  any  deep  knowledge  or  skill  in  APL 
programming  and  computer  operations.  On  the  other  hand, 
the  sophistication  of  the  advanced  user  is  taken  into 


consideration. 
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Further  use  of  the  package  will  show  if  these 
goals  have  been  achieved  successfully.  The  first  attempt 
cannot  be  the  final  and  optimal  solution.  Reprogramming 
or  writing  additional  APL  functions  by  the  maintenance 
programmer  may  be  needed.  Therefore,  only  the  basic 
functions  are  included  in  the  package  and  other  desirable 
ones  are  suggested  as  future  work  or  as  possible 
improvements . 

The  reader  interested  primarily  in  practical  use 
of  the  package  should  start  with  Section  4.10  and  proceed 
to  Chapters  3,  6  and  others.  Section  4.10  includes  a 
detailed  scenario  illustrating  use  of  the  package,  with 
many  examples  of  the  outputed  maps. 
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2.  STATE  OF  THE  ART  IN  COMPUTER  CARTOGRAPHY 

2 . 1  General  Overview 

Computer  mapping  is  a  part  of  the  larger  field  of 
computer  cartography.  It  consists  of  using  computers  to 
produce,  manipulate,  generate,  combine  and  perform  logical 
and  arithmetic  operations  on  spatially  distributed  data 
represented  by  maps.  "It  is  probably  the  most  important 
advance  made  in  the  field  of  cartography  over  the 

past  twenty-five  years  and  the  potential  for  the  future  is 

v 

enormous."  (18.) 

However,  computer  cartography  is  still  passing 
through  a  stage  of  infancy.  Important  general  aspects  to 
be  considered  are:  (a)  data  and  associated  problems, 

(b)  hardware,  (c)  software,  (d)  methods  and  theoretical 
justification,  and  (e)  applications  and  people  working 
in  the  field. 

2.1.1  Data 

Data  used  in  computer  cartography  are  very 
voluminous  and  collected  on  a  wide  basis  for  various 
purposes.  Many  of  the  existing  data  sources  were  estab¬ 
lished  a  long  time  ago,  and  data  was  obtained  using  various 
methods  ranging  from  aerial  photography  to  entry  by  hand 
tracing  of  an  existing  picture. 
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The  computer’s  most  important  contributions  to 
the  field  are  automatic  digitalization  and  the  organization 
of  large  collections  of  data  into  "data  banks".  An 
example  of  the  further  development  of  automated  data 
collection  is  the  use  of  satellite  photographs  of  the 
earth’s  surface. 

The  first  symposium  on  geographical  information 
systems  was  held  in  Ottawa  in  September  1970.  It  dealt 
with  the  acceptance,  storage,  manipulation,  retrieval  and 
display  of  geographical  and  environmental  data.  Attempts 
were  made  to  summarize,  standardize,  and  direct  the  disjoint 
efforts  in  this  field. 

2.1.2  Hardware 

Hardware  used  in  computer  cartography  is  still 
very  inefficient  for  that  purpose.  The  great  majority 
of  computer  maps  have  been  produced  on  the  line  printer 
because  of  its  ready  availability,  speed,  and  low  cost. 

Since  the  line  printer  was  not  designed  to  produce  maps, 
many  conveniences  have  to  be  eliminated.  For  example, 
continuous  lines  cannot  be  drawn,  and  printed  characters 
must  be  used. 

Mechanical  plotters  are  more  appropriate  devices 
for  map  making.  Flatbed  plotters  and  drum  plotters  are 
used.  Flatbed  plotters  offer  a  large  plotting  area  but 
are  relatively  slow.  Multicoloured  plots  can  be  produced. 
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Accuracy  in  most  cases  is  satisfactory.  Drum  plotters 
move  the  pen  in  one  direction  and  rotate  the  paper 
perpendicular  to  the  pen’s  moving  direction.  Plots  can 
therefore  have  infinite  dimension  in  one  direction. 

Slowness  of  the  plotters  is  still  the  most  limiting  factor. 

Electronic  plotters  use  a  cathode  ray  tube  for 
image  generation,  and  then  photograph  it  to  produce  hard 
copy  output.  Maps  can  be  produced  rapidly  and  accurately, 
but  the  surface  area  is  small,  and  the  CRT  is  still  not 
available  on  a  wide  basis.  Other  relatively  new  devices 
are  available,  for  example,  a  plotter  using  a  xerox- copying 
principle,  but  are  not  yet  widely  used. 

2.1.3  Software 

Computer  cartography  software  is  still  limited  to 
a  few  techniques  and  widely  accepted  systems,  but  progress 
is  expected  soon.  This  is  the  part  of  the  field  that  most 
people  are  working  on.  Some  of  the  developed  programs  and 
systems  will  be  discussed  later  in  this  section. 

2.1.4  Methods 

The  worst  situation  is  in  theoretical  methods 
development  and  justification.  Very  few  people  are  working 
in  this  area.  The  majority  of  theoretical  computer  carto¬ 
graphers  are  concerned  with  the  application  of  statistical 
methods.  The  further  development  of  quantitative  computer 
cartography  and  theoretical  geography  will  hopefully  facilitate 
the  development  of  theoretical  approaches  and  methods. 
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2.1.5  People 

People  working  in  computer  cartography  come  from 
different  fields  of  application.  The  mostly  develop  programs 
for  the  meaningful  representation  in  map  form  of  their 
field’s  data.  Yet  some  general  systems  have  been  developed, 
and  centers  for  computer  cartography  have  been  established  and 
have  gained  wide  acceptance  and  reputation  in  the  field. 

In  1966,  H.  T.  Fisher  established  the  Laboratory 
for  Computer  Graphics  and  Spatial  Analysis  at  Harvard 
University.  Since  then  the  laboratory  has  produced  many 
programs  for  computer  mapping  that  are  applicable  to  many 
fields.  Another  center  is  at  the  University  of  Michigan  in 
the  Department  of  Geography,  headed  by  W.  R.  Tobbler. 

J.  L.  Pfalz,  now  at  the  University  of  North  Carolina,  and 
A.  Rosenfield  at  the  University  of  Maryland  have  done  signi¬ 
ficant  work.  In  Canada,  D.  R.  F.  Tomlinson  and  D.  H.  Douglas 
at  Carleton  University  in  Ottawa,  A.  R.  Boyle  at  the  University 
of  Saskatchewan,  and  T.  Peucker  at  Simon  Fraser  University 
in  Burnaby,  B.C.,  have  been  active  in  the  field.  Bengtsson 
and  Nordbeck  at  the  University  of  Lund,  Sweden  are  currently 
contributing  to  the  field.  Experimental  work  is  being 
carried  on  in  the  Universities  of  Edinburgh  and  Glasgow. 

Some  of  the  works  done  will  be  briefly  discussed 
later  in  this  thesis,  as  well  as  some  of  the  many  applications 
that  can  be  found  in  literature. 


-8- 


2 . 2  Theoretical  Work,  Formal  Definitions 

In  computer  cartography,  very  little  has  been 
done  on  theoretical  approaches  to  the  field.  This  is 
probably  because  few  people  have  been  working  exclusively 
on  computer  cartography,  and  also  because  they  are  more 
interested  in  developing  useful  programming  systems  for 
real-life  applications. 

2.2.1  Contour  Lines 

Much  work  has  been  done  on  contour  line  representa¬ 
tions  of  maps.  Individuals  active  in  this  area  include 
S.  F.  Morse,  H.  Freeman,  and  B.  W.  Boehm. 

Many  real  life  functional  relationships  are  so 
irregular  that  well  defined  mathematical  functions  are  not 
suitable  for  their  representation.  Consider  the  class  of 
functions  z=f(x,y)  .  An  example  of  such  a  function  is  the 
altitude  of  a  given  area  of  terrain.  Other  examples  can 
be  found  in  meteorology,  oceanography,  demography,  social 
sciences  and  engineering.  The  function  f  can  only  be 
defined  by  listing  the  set  of  doublets  ((x,y),z)  in 
tabular  form.  Instead  of  representing  the  data  in  tables, 
curves  of  equal  values  of  the  dependent  variable  z  ,  called 
isarithms,  can  be  drawn  in  two  dimensions,  the  dimensions 
representing  the  independent  variables  x  and  y  .  These 
contour  lines  make  up  a  contour  map.  This  method  is  of 
great  help  to  the  human  eye  since  certain  visual  patterns 
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o£  the  map  indicate  properties  of  the  terrain.  For  example, 
the  density  of  the  contour  lines  gives  an  indication  of  the 
magnitude  of  the  slope.  Morse  S.  P.  (1.)  defines  contour 
lines  as  follows:  Let  the  altitude  function  z=f(x,y)  be 
continuous  and  single-valued,  that  is,  there  are  no  vertical 
or  overhanging  cliffs.  Assign  elevation  values  e^  ,  such 
that  -  e^  =  E  where  E  is  a  constant  altitude 

increment  and  z  <e. <z  for  all  i  =  l,2,...,n  . 

For  each  i  there  exists  a  contour  line  of  elevation  e^ 
that  separates  the  points  of  equal  or  higher  altitude  from 
the  points  of  equal  or  lower  altitude.  Out  of  these  formal 
definitions  come  useful  concepts  of  the  positive  and  negative 
contour  lines  at  the  edges  of  plateaus.  Further  definitions 
are  the  normal  and  degenerate  contour  lines  (both  positive 
and  negative),  or  maximum  and  minimum  contour  lines. 

Using  the  contour  line  definition  various 
properties  regarding  elevation  types  of  the  points  on  or 
near  the  contour  lines  are  derived.  Directions  are 
assigned  to  the  lines  such  that  their  high  sets  are  always 
on  the  left. 

Two  contour  lines  are  adjacent  if  they  can  be 
connected  by  a  line  that  does  not  cross  another  line.  From 
this,  various  types  of  adjacency  sets  are  developed  which 
can  be  used  to  generate  a  graph  of  the  contour  map. 

A  supplementary  definition  of  a  contour  line 
given  by  Boehm  (2.)  is  that  g(Xj,y^)  approximates 
z(x,y)  over  the  region  R  where  ^xj  »Tp)  in  R  are 
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chosen  such  that  g(x^,yk)  =  e±  .  The  value  at  any  other 

point  is  found  either  by  taking  the  nearest  (x^ ,yk)  and 

its  value  g(x.,yi)  ,  or  by  use  of  interpolation.  If  the 

J  k 

(x-,yk)  are  discrete  points  one  obtains  a  grid  representation. 

Some  more  definitions  are  given  by  Morse.  It 
is  assumed  that  the  map  can  be  divided  into  finite 
neighbourhoods.  These  then  are  used  to  define  elevation 
types  which  classify  subcritical  neighbourhoods  (the  type 
does  not  change  inside  them)  according  to  the  elevation 
of  their  points.  Any  combination  of  elevation  equal  to, 
greater  than  or  smaller  than  the  given  point  make  up  seven 
elevation  types.  The  latter  are  the  major  tools  in  the 
proofs  given  for  the  various  properties  of  contour  lines. 

Closed  contour  lines  may  form  concave  shapes 
resulting  in  caves,  cavoids  and  doors. 

A  simple  maze  in  a  contour  map  is  defined  as 
the  focus  of  all  straight  lines  that  intersect  two  disjoint 
closed  contour  lines. 

2.2.2  Pictures  and  Picture  Processing 

The  work  done  by  A.  Rosenfeld  at  the  University 
of  Maryland  does  not  deal  particularly  with  maps  as  such, 
but  with  more  generally  defined  objects  called  "pictures’' 
or  "picture  functions".  The  definition,  however,  includes 
maps.  Many  of  the  results  of  that  work  could  be  used  in 
theoretical  work  in  computer  cartography.  Pictures  are 
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defined  in  that  work  as  mathematical  objects,  and  the  use 
of  mathematical  operations  on  pictures  is  justified  by 
showing  that  picture  functions  are  "indistinguishable" 
from  functions  which  are  well  behaved  in  various  senses. 

Types  of  picture  processing  problems  are  outlined 
and  considered.  Areas  covered  include  the  following: 

a)  the  picture  as  an  information  source  and  their 
efficient  encoding. 

b)  the  approximation  of  pictures  -  sampling  and 
quantization  techniques. 

c)  position  variant  operations  on  pictures  and 
their  implementation. 

d)  picture  properties  useful  for  pictorial  pattern 
recognition. 

e)  "figure  extraction",  figure  properties. 

f)  picture  description  and  "picture  languages". 

In  the  other  work  (4.), proofs  for  border  following  algorithms 
of  simply-connected  objects  are  given. 


2 . 3  Map  Representation  Techniques 

To  represent  a  map  in  the  computer,  various 
methods  have  been  successfully  used  and  they  are  briefly 
discussed  in  the  following.  No  single  representation  is 
optimal  for  all  applications.  Criteria  for  the  comparison 
in  terms  of  the  tradeoff  between  storage  space  and 
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computation  speed  have  been  suggested  (2.). 

2.3.1  Contour  Lines  and  Chain  Encoding 

Contour  lines  can  be  used  to  represent  the  necessary 
information.  This  is  a  particularly  economical  way  to 
represent  a  map  if  the  chain  encoding  method  developed  by 
H.  Freeman  is  used  (5., 6.).  The  values  for  f  are  no 
longer  given  for  each  point  of  the  terrain,  but  only  for 
a  few  points  selected  in  a  specific  way. 

Chain  encoding  is  a  general  technique  applicable 
also  to  the  encoding  of  any  irregular  curve  that  cannot 
be  given  conveniently  as  a  mathematical  function.  To 
describe  the  chain  encoding  scheme  briefly,  consider  the 
curve  shown  in  Figure  1.  A  fine  square  grid  is  assumed 
to  be  superimposed  on  the  curve.  The  curve  is  traced  and 
intersections  between  the  curve  and  the  grid  are  noted. 

At  each  intersection,  the  distances  to  the  two  nodes 
surrounding  this  intersection  are  compared.  The  closer 
of  the  two  nodes  is  marked  as  a  chain  point  and  becomes  a 
part  of  the  straight  line  approximation  used  to  represent 
the  curve.  Thus  X,  Z,  V,  T  are  successive  chain  points. 
Obviously  the  X^  +  -^  chain  point  can  lie  only  at  one 

of  the  eight  nodes  surrounding  the  X^th  chain  point. 

This  is  suggested  by  the  simple  encoding  scheme  illustrated 
in  Figure  2  where  eight  possible  directions  are  numbered  in 
a  clockwise  fashion  starting  from  the  right  horizontal. 
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Figure  1 


Figure  2 
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The  chain  of  the  given  curve  is  thus  formed  by 
sequences  of  octal  numbers.  The  curve  represented  in 
Figure  2  has  been  encoded  into  the  chain  001102243  . 
Depending  on  the  grid,  approximation  can  be  made  accurate 
enough  for  an  intended  purpose.  The  chain  encoding  repres¬ 
entation  of  a  contour  line  permits  the  development  of  simple 
algorithms  for  determination  of  chain  length,  determination 
of  maximum  excursions  along  the  x  and  y  axes  (called 
width  and  height  of  the  chain),  determination  of  the  area, 
determination  of  residue  of  a  chain,  etc.  The  synthesis 
and  manipulation  of  chain  encoded  curves  can  also  be 
performed  without  major  difficulties.  Synthesis  is  defined 
as  the  generation  of  a  line  drawing  to  satisfy  certain 
specifications.  Manipulation  is  the  transformation  of  a 
line  drawing  in  a  specific  manner  into  another  line  drawing. 

2.3.2  Polygon  Representation 

Another  method  of  representation  is  the  digitized 
map  where  the  map  is  represented  in  the  form  of  contour 
lines  approximated  by  polygons.  This  representation  can  be 
advantageous  if  areas  and  longitudes  of  contour  lines  are 
frequently  needed  or  calculated. 

2.3.3  Graph  Representation 

Contour  maps  can  be  represented  as  a  graph  (2.). 
The  closed  curves  become  links  and  the  areas  between  become 
nodes.  The  graph  is  directed  from  the  external  contour 
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lines  toward  the  internal  ones  (5.)*  Any  contour  line 
can  only  be  internal  to  one  other  contour  line;  hence, 
only  one  link  can  enter  a  node.  The  graph  therefore  is 
a  tree,  the  root  corresponding  to  the  map  boundary. 

2.3.4  Functional  Representation 

A  representation  of  a  map  using  mathematical 
functions  has  been  worked  out  by  McMenamin  (7.).  A 
particular  application  describes  a  simulation  of  views 
as  seen  by  astronauts  in  a  Lunar  Surface  Roving  Vehicle. 

The  terrain  is  represented  by  a  series  of  Lagrange  poly¬ 
nomials.  Only  the  coefficients  of  the  polynomials  for  each 
grid  square  are  stored.  Some  of  them  are  shared  between 
two  or  more  grid  squares.  The  coefficients  represent 
partial  derivatives  of  the  altitude  function  with  respect 
to  the  coordinates  up  to  the  4th  order. 

2.3.5  Grid  Point  Representation 

A  map  can  be  represented  at  grid  points  by  the 
values  of  the  function  f(x,y)  .  These  may  be  spaced 
uniformly  or  their  spacing  may  be  variable  depending  on 
the  irregularity  of  the  terrain.  The  visual  effect  of  such 
a  map  is  not  as  graphic  as  the  contour  lines.  This 
representation  is  used  in  civil  engineering. 

The  grid  method  with  a  differential  representation 
can  be  used.  Only  variations  from  point  to  point  are  stored 
with  the  absolute  value  of  one  reference  point  added. 
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Maps  can  be  also  represented  with  only  coarse 
grid  points  as  reported  in  the  work  of  Boehm  and  Rieber 
(8.)*  The  values  for  points  in  between  may  be  generated 
by  interpolating  from  the  neighbouring  points.  Instead 
of  using  a  linear  interpolation  the  probability  density 
function  (PDF)  of  variations  for  the  particular  terrain 
is  used.  Distinct  PDF’s  are  used  for  the  X  and  Y 
direction.  The  use  of  joint  instead  of  marginal  PDF's 
is  suggested  for  further  work  (2.). 

In  this  thesis,  maps  are  synthesized  from  a  set 
of  independent  information  units  concerning  an  area,  called 
base  maps.  A  base  map  is  a  type  of  grid  map.  Base  maps 
can  be  combined  logically  or  mathematically.  The  method 
will  subsequently  be  discussed  in  detail. 

2 . 4  Programs  and  Systems:  Developed  and  Used 

Many  developed  systems  and  programs  have  been 
reported.  Some  originating  at  universities  have  found 
wider  application  because  of  their  generality  than  others 
developed  for  special  purposes  at  companies,  agencies, 
and  governmental  organizations. 

Only  a  few  of  them  will  be  briefly  mentioned  here. 
General  information  about  most  existing  systems  can  be 
found  in  the  literature  (9.). 
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2.4.1  SYMAP 

The  best  known  and  probably  most  widely  used 
computer  mapping  package,  SYMAP,  was  developed  in  1963  by 
H.  Fisher  at  Northwestern  University.  Fisher  later 
established  the  Laboratory  for  Computer  Graphics  and 
Spatial  Analyses  at  Harvard  University.  The  laboratory 
devotes  much  of  its  work  to  the  maintenance  and  improvement 
of  SYMAP  and  to  the  development  of  other  programs  used  in 
computer  mapping. 

The  latest  version  of  SYMAP  offers  three  basic 
formats  for  output  with  some  35  options  available.  Up  to 
ten  levels  of  shading  can  be  obtained  by  using  the  over¬ 
printing  technique.  The  formats  are:  contour  (isoline) 
map,  conformant  (chloropleth)  map,  and  proximal  map. 

The  contour  (or  isoline)  map  consists  of  closed 
curves  connecting  the  points  with  the  same  numeric  value. 

The  contour  map  is  suitable  for  representation  of  continuous 
information. 

Conformant  (or  cloropleth)  maps  are  based  on  the 
conformance  to  the  boundaries  of  the  data  zone.  This  type 
is  suitable  for  data  whose  areal  limits  are  of  significance 
and  whose  representation  as  a  continuous  surface  is 
inappropriate . 
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Proximal  maps  are  based  on  proximity  to  a  data 
point.  Each  character  location  on  the  output  map  is 
assigned  a  value  to  the  nearest  data  point.  This  type  of 
map  is  particularly  useful  when  the  data  is  of  a  qualitative 
rather  than  a  quantitative  nature  and  when  the  categorizing 
of  data  cannot  be  assumed. 

Two  examples  of  SYMAP  output  are  shown  on  the 
Figures  3.  and  4.  SYMAP  output  can  show  only  overall 
patterns.  The  maps  are  not  a  precise  presentation  of  the 
data.  The  output  also  cannot  be  readily  interpreted  by 
persons  unfamiliar  with  SYMAP.  Nonetheless,  SYMAP  is  very 
widely  used. 

There  are  many  mapping  programs  with  output 
similar  to  SYMAP  written  to  satisfy  specific  requirements 
that  cannot  be  done  by  SYMAP.  These  programs,  also 
developed  at  the  Laboratory  for  Computer  Graphics  and 
Spatial  Analyses  at  Harvard,  are:  RUMOR,  GRID,  OTOTROL, 
LOKAT,  GRASP,  OBLIX  1,  and  SYMVU. 

GRID  is  used  to  map  large  quantitites  of  data 
collected  on  the  basis  of  a  rectangular  coordinate  grid. 

OTOTROL,  LOKAT,  and  GRASP  are  programs  for  three 
dimensional  applications,  and  can  be  applied  effectively 
in  architecture. 

OBLIX  1  produces  maps  of  surfaces  with  contour 
lines  included  and  provide  indications  of  background 
texture  from  an  input  matrix. 


Fig. 3. 
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SYMAP’s  wide  use  is  due  to  its  relatively  early 
appearance,  good  advertising,  good  documentation,  and  ease 
in  using,  and  also  because  of  its  ability  to  produce  maps 
on  a  line  printer  in  an  inexpensive  way. 

2.4.2  Other  Systems 

Similar  programs  to  SYMAP  were  also  developed 
by  Bengtsson  and  Nordbeck  at  the  University  of  Lund,  Sweden, 
and  also  by  G.  M.  Gaits  in  Great  Britain  (LINMAP  and  COLMAP) . 

At  the  University  of  Maryland,  the  Map  Analyses 
System  (MANS)  has  been  developed  by  J.  L.  Pfalz  (10.).  For 
the  purpose  of  the  MANS  system  a  map  is  defined  to  be  a 
digital  picture  in  which: 

a)  each  connected  component  has  an  individual 
significance  for  the  analyses  of  the  picture, 

b)  the  number  of  connected  components  is  "small" 
relative  to  size  of  the  picture. 

(Typically  it  might  be  expected  that  an  n  x  n  map  contains 
fewer  that  2n  distinct  connected  components) .  Conventional 
maps  and  line  drawings  will  usually  satisfy  these  require¬ 
ments.  The  analyses  of  a  picture  involves  non-pictorial 
information  in  the  form  of  statements  about  the  nature  of 
objects  which  appear  in  the  picture.  In  practice  these 
statements  may  not  be  given  explicitly  but  in  the  form  of 
symbols.  For  example  +  implies  "This  is  a  church".  The 
analyses  of  a  picture  or  map  includes: 
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a) 

classification 

of 

the 

entire 

picture , 

b) 

isolation  and 

clas 

s if ication 

of  specific 

of  a  picture, 

c) 

description  of 

geometrical  properties  of 

individual  sub 

sets 

of 

the  picture. 

The  output  of  the  MANS  system  is  in  the  form  of  displayed 
pictures  created  from  the  original  picture  together  with 
statements  containing  information  about  the  displayed 
picture.  The  display  in  the  present  implementation  is 
drawn  on  the  plotter  but  should  ideally  occur  on  a  CRT. 

In  the  examples  shown  in  Figures  5.  and  6.,  MANS 
has  found  specified  entities  and  shaded  them. 

Other  systems  and  programs  are  now  mentioned  very 

briefly : 

The  Canada  Geographic  Information  System  is 
utilized  in  planning  land  resources  use. 

The  U.S.  Geological  Survey  System,  developed  at 
the  U.S.  Geological  Survey  in  Washington,  is  intended  to 
provide  base  mapping  for  the  United  States,  and  will  be 
completed  in  10  years.  Hardware  and  Software  are  presently 
being  developed  to  merge  new  map  information  onto  the  base 
maps . 

The  Remote  Sensing  and  Environment  System  has 
been  developed  at  Cornell  University.  The  basic  approach 
for  this  system  is  to  obtain  the  map  information  by  air 
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Figure  5. 
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Figure  6. 
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photography  or  by  a  field  survey,  to  manually  record, 
keypunch  and  store  the  data.  Twenty-five  generalized 
programs  are  used  for  retrieval,  and  give  tabular  or 
graphic  output  on  a  printer.  The  system  is  based  on  grid 
manipulation. 

The  Munich  Air  Pollution  Survey  System,  developed 
at  the  University  of  Wurzburg,  employs  remote  centres 
feeding  a  data  bank.  From  this  a  meteorological  forecast 
is  made.  The  programs  for  plotting  forecast  maps  are 
still  in  the  stage  of  development. 

The  AUTOMAP  -  Automatic  Mapping  System  developed 
at  the  C.I.A.  in  Washington  has  been  operational  since  1967, 
and  can  produce  any  form  of  line  or  point  data. 

2 . 5  Interactive  Mapping 

Although  interactive  computer  mapping  has  been 
mentioned  in  several  articles,  there  have  been  very  few 
attempts  to  develop  interactive  mapping  systems. 

In  cooperation  with  the  Canadian  Hydrographic 
Services,  the  Electrical  Engineering  Department  of  the 
University  of  Saskatchewan  has  been  examining  the  feasib¬ 
ility  of  various  systems  and  methods.  As  reported  by 
A.  R.  Boyle  at  a  seminar  sponsored  by  the  Data  Systems 
Section  of  the  Radio  and  E.  E.  Division  of  the  National 
Research  Council  (May  31,  and  June  1,  1971),  this  work 
involves  hardware  and  software  development.  It  is 
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especially  concerned  with  detailed,  high  accuracy  mapping 
called  Computer  Aided  Compilation  (CAC) .  The  system  is 
designed  to  allow  off-line  preparation  of  a  display  tape 
in  a  large  computer,  followed  by  interactive  decision-making 
carried  out  at  a  display  driven  by  a  small  computer  fitted 
with  a  low  cost  disc.  The  Tektronix  611  storage  display  is 
fitted  with  both  a  keyboard  and  a  type  of  ’’pointer”,  a 
"non-stored”  spot  which  is  moved  over  the  screen  by  means 
of  a  joystick,  rolling  ball,  ’’mouse"  or  input  tablet. 

A  typical  use  of  the  system  is  as  follows: 
Hydrographic  data  on  coastlines  sandings,  and  water  contours 
is  entered  into  the  large  computer  and  a  tape  is  prepared. 

The  user  can  then  interact  with  a  map  representing  this  data. 
He  can  add  or  erase  graphic  features,  change  their  method 
of  representation,  translate  or  scale  the  picture,  and  add 
new  topographic  data.  When  work  on  any  area  has  been 
completed,  an  output  tape  is  generated,  containing  the 
modification  data  as  well  as  original  data.  This  is 
returned  to  the  off-line  computer  and  modifications  are 
incorporated . 

Another  interactive  system,  already  described, 

is  MANS. 

In  this  thesis  an  interactive  approach  is  also 
attempted.  To  make  a  clear  distinction  between  the  use 
of  the  full-scale  CRT  interactive  system  and  this  mode  of 
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interaction  using  a  remote  APL  terminal,  the  latter  will 
be  called  conversational. 

2 . 6  Applications 

Computer  mapping  has  been  widely  used  in  various 
fields  and  is  a  truly  multidisciplinary  tool.  Any  type  of 
data  which  has  spatial  significance  can  be  illustrated, 
compared,  aggregated,  manipulated,  and  analysed  by  a  variety 
of  programs  producing,  as  a  final  output,  information  in  the 
form  of  a  map. 

Obviously  cartography  itself  is  (or  should  be) 
the  principal  field  of  computer  application,  especially 
where  accuracy  and  speed  are  of  importance  and  the 
complexity  of  the  map  relatively  small.  Future  prospects 
are  positive.  The  demand  for  automated  mapping  is  increasing 
with  efforts  to  gather  data  for  geographical,  urban, 
geological  and  other  data  banks. 

In  landscape  architecture,  city  and  regional 
planning,  urban  design  and  other  environmental  sciences 
there  is  an  increasing  demand  to  employ  computer  oriented 
techniques.  It  can  facilitate  visual  understanding,  and 
thus  open  wider  possibilities  for  planners,  researchers 
and  designers. 

A  few  of  the  applications  will  now  be 


mentioned. 


-27- 


A  project  which  is  summarized  as  a  computer  map 
atlas  illustrates  the  application  of  SYMAP  in  environmental 
sciences.  MA  Computer  Atlas  of  Ottawa-Hull"  (18.)  was 
elaborated  by  D.  R.  F.  Taylor  and  D.  H.  Douglas  at  the 
University  of  Ottawa  on  an  IBM  360/65  Computer,  using  also 
Autotrol  Digitizer  Model  3700  as  the  main  graphic  input 
tool.  The  maps  assembled  dealt  with  age  structure,  ethnic 
factors,  political  responses,  and  economic  distributions. 
Figures  3.  and  4.  are  taken  from  this  Atlas. 

The  paper  of  S.  H.  Brownstein  (12.)  describes  two 
projects  established  to  develop  and  demonstrate  computer 
based  techniques  for  creating  and  maintaining  a  geograph¬ 
ically  oriented  urban  data  base,  and  retrieving,  analysing 
and  displaying  meaningful  information.  Methods  have  been 
developed  to  facilitate  specifying  geographic  areas  of 
interest  and  desired  information,  and  then  to  retrieve 
relevant  records. 

The  Urban  Planning  Data  Systems  Project  at  the 
University  of  Washington  deals  also  with  similar  problems. 

M.  Y.  Kevany  describes  selected  concepts  and 
techniques  of  data  processing  system  in  the  document 
"An  Information  System  for  Urban  Transportation  Planning" 
(19.)  and  in  "BATSC  Approach  Technical  Memorandum  for  System 
Development  Corporation".  The  system  is  in  use  at  Bay  Area 
Transportation . 
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The  range  o£  other  applications  is  very  wide  and 
only  a  few  interesting  applications  will  be  mentioned  here. 

Contour  map  analysis  has  been  used  in  finding  the 
flight  path  for  a  given  terrain  profile  and  a  direction  of 
flight  of  an  aircraft  (5.).  Terrain  has  been  simulated  for 
lunar  vehicles  and  low  flying  airplanes  (7.). 

Transitive  maps  have  been  used  in  demographic 
studies  by  Nordbeck  and  Rystedt  (13.).  The  Range  Map 
Function  g(x,y)  is  defined  as  the  radius  of  the  circle 
centered  at  (x,y)  ,  that  is  required  to  enclose  a  standard 

size  population.  The  Carriage  Map  function  f(x,y)  is 
the  cost  of  transportation  of  a  standard  amount  of  goods 
to  a  standard  size  population  that  is  closest  to  (x,y)  . 

Equal  values  of  g  and  f  are  plotted  and  generate  a 
contour  map. 

An  interesting  attempt  has  been  made  by  B.  Cornwall 
and  D.  H.  Robinson  (14.).  The  paper  describes  their  work 
in  animated  mapping  with  computers,  using  the  CRT  with 
fixed  intensity  beam  as  the  drawing  "tool".  The  image  was 
recorded  by  a  shutter-less  cine  camera.  Drawing  of  maps 
by  computer  with  slight  successive  changes  results  in  an 
animated  map.  Examples  are  given  from  several  fields. 

There  are  many  other  applications  in  various 
fields,  for  example:  physics  (X-ray  crys talography) , 
oceanography,  economics  (market  analysis  study),  highway 
design  (determination  of  an  optimal  highway  route  through 
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mountainous  terrain) ,  geophysics  (mapping  of  a  geomagnetic 
field) ,  geology  (surface  determination  for  mineral 
exploration  using  contour  maps) ,  sociology  (population 
changes  and  population  density  maps),  statistics,  forestry, 
traffic,  civil  engineering  (determination  of  terrain  to 
be  flooded  by  a  dam  placed  at  a  particular  location,  terrain 
visibility  from  a  particular  top  of  the  mountain  for  TV 
Tower  locations),  military  applications  and  many  others. 

Two  illustrations  of  the  use  of  maps  are  given  in 
Figures  7.  and  8.  Figure  7.  shows  the  Boston  industry 
and  traffic  flow.  Figure  8.  shows  the  population  Graph- 
Map  of  the  U.S.A.  Both  are  produced  by  the  Laboratory 
for  Computer  Graphics  and  Spatial  Analysis  at  Harvard 
University. 
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Figure  7. 


Figure  8 
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3.  APL  BASED  PACKAGE  DESIGN  CONCEPTS 
3 . 1  Base  Maps 

To  represent  spatially  distributed  environmental 
factors,  a  binary  matrix  representation  is  used.  Each 
element,  X^  ,  of  the  matrix  corresponds  to  the  grid 
point  with  coordinates  (i,j)  .  If  the  desired  factor  is 
present  at  (i,j)  on  the  grid,  the  matrix  element 
is  set  to  1  .  Otherwise,  it  is  set  to  0  .  This  matrix 
representation  of  certain  characteristics  of  an  area  is 
called  a  base  map.  Thus,  for  each  factor  employed,  one 
base  map  is  generated. 

The  input  of  base  maps  will  be  discussed  in  the 
following  chapter.  It  is  not  necessary  to  input  all  the 
grid  points  of  the  area  containing  the  desired  factor. 

Like  the  chain  encoding  technique,  only  a  set  of  points 
which  determine  the  curve  that  encloses  the  area  are 
required.  If  there  is  no  area,  i.e.,  only  discrete  points 
are  used,  all  points  must  be  input  when  generating  the 
base  map.  Discrete  points,  however,  are  not  usually  found 
when  dealing  with  environmental  factors  (land,  lakes, 
rivers,  etc.).  These  are  almost  always  contained  within 
an  area. 

There  are  certainly  cases  when  binary 
representation  does  not  satisfy  the  needs  of  the  user,  as, 
for  example,  in  the  representation  of  the  cost  of  land  parcels. 
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All  cost  values  should  be  represented  in  one  map.  The 
functions  in  the  package  are  not  restricted  to  the  handling 
of  binary  maps,  but  APL  storage  (the  libraries  and  workspace) 
imposes  serious  limitations  in  these  cases. 

Clear  workspace  has  31868  bytes.  Integers  take 
4  bytes,  reals  take  8,  literals,  1,  and  logical  or  binary 
data,  1  bit.  Thus,  the  real  value  matrix  with  dimensions 
50  x  70  used  in  this  thesis  would  use  up  more  than  2/3 
of  the  available  workspace.  In  order  to  avoid  such  a 
situation,  each  value  can  be  represented  on  a  different 
base  map.  More  detailed  discussion  of  this  problem 
will  occur  in  subsequent  chapters. 

Another  problem  is  the  size  of  the  base  maps. 

If  too  large,  base  maps  would  result  in  WS  FULL  messages 
whenever  they  were  employed  in  an  expression,  or  used  with 
certain  functions.  Tricky  manouevring  of  the  developed 
functions  could  sometimes  avoid  this  difficulty,  but  the 
overhead  of  such  manipulations  could  be  considerable.  The 
50  x  70  base  map  was  chosen,  since  it  is  sufficiently  big 
to  represent  an  area  and  fits  the  standard  form  used  for 
xeroxing.  Thus,  reproduction  can  easily  be  made  from  the 
original  obtained  on  the  terminal. 

This  representation  certainly  has  both  advantages 
and  disadvantages.  Its  main  advantage  is  its  simplicity. 
Also,  because  of  the  binary  representation,  all  logical 
operations  that  APL  provides  are  immediately  applicable. 
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Super  imposing  environmental  factors  can  be  achieved  with 
logical  operations  on  base  maps.  Arithmetic  operations  are 
also  feasible  since  the  base  map  contains  logical  (0,1)  , 

algebraic  (0,1  in  decimal  system),  and  binary  (0,1  in 
binary  system)  values  at  the  same  time,  depending  only  on 
the  semantics  of  the  APL  expression. 


Insertion  or  deletion  of  points  is  done  in  a 
straightforward  way,  using  the  APL  assignment  operator. 
Areas  too  can  be  changed  using  the  same  procedure  as  at 
input . 

Input  of  areas  can  be  done  in  the  simplest  way 
(similar  to  the  chain  encoding  technique)  by  entering  only 
the  vector  of  points  that  define  the  contour  of  the  area. 


The  main  disadvantage  of  the  base  map  representation 
is  that  it  takes  up  too  much  room.  The  workspace  is  crucial. 
The  actual  size  used  is  actually  a  compromise.  Other 
representations  would  require  some  conversions  that  would 
take  considerable  time  whenever  logical  operations  are  to 
be  performed.  For  example,  a  chain  encoded  representation 
using  integers  would  not  save  as  much  in  storage  costs  as 
it  would  waste  in  the  cost  of  computing  time. 


APL  PLUS,  an  extension  of  APL,  allows  the  user 
to  organize  his  own  files  using  very  simple  APL- like 
instructions,  and  to  store  the  base  maps  in  the  file. 
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In  this  work,  the  current  version  of  the  APL  PLUS 
File  Subsystem,  developed  by  I.  P.  Sharp  Associates  Limited, 
Toronto,  was  used.  This  version  is  commercially  available. 
It  gives  to  APL/ 360  two  fundamental  new  capabilities: 


1)  the  entry,  retention,  and 
control)  of  much  more  data  than  can 

2)  the  storing  and  selective 
base  by  a  group  of  users. 


access  (under  programme 
be  stored  in  an  APL  workspace 
modification  of  a  data 


A  file  in  this  system  is  a  collection  of  data  organized  into 
components.  A  component  of  a  file,  like  a  variable  in  a 
workspace  may  be  any  APL  value - s calar ,  vector,  matrix,  or 
higher-dimensional  array.  The  only  limit  on  the  size  of  a 
component  is  that  it  must  fit  in  a  workspace.  Obviously, 
this  concept  of  a  file  system  is  very  convenient  for  base 
map  storage.  For  each  entity  shown  on  a  set  of  base  maps, 
a  file  can  be  organized  with  base  maps  as  its  components. 

The  relative  overhead  of  APL  file  information  decreases  as 
the  component  becomes  longer. 


The  file-sharing  concept  opens  up  a  wide  range 
of  possibilities  for  more  than  one  user.  A  file  owner  may 
specify  access  control  for  a  file  to  designate  what  other 
persons  (account  numbers)  may  access  the  file  and  which 
file  operations  they  may  perform. 
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3 . 2  Value  Maps 

The  matrices  whose  elements  are  values  (integer 
or  real)  are  allowed  in  this  system.  They  are  called 
value  maps.  Value  maps  may  result  from  arithmetic  operations 
on  base  maps,  from  the  user's  definition  or  input,  or  from 
retrieval  from  a  file  in  which  they  have  been  stored. 

One  way  in  which  value  maps  can  be  used  results 
essentially  in  their  conversion  to  binary  maps.  For  example, 
if  the  value  map  named  MAPI  is  used,  and  points  with  value 
A  are  of  particular  interest,  the  expression  (MAP1=A) 
can  be  used.  This  expression  results  in  a  matrix  that  has 
the  dimensions  of  MAPI  ,  with  a  value  of  0  for  all 
elements  unequal  to  A  ,  and  a  value  of  1  for  all 
equivalent  elements.  This  is  nothing  more  than  a  binary  map, 
i.e.,  a  base  map  showing  all  points  that  have  value  A  . 

It  should  definitely  be  possible  to  output  the 
value  maps;  therefore,  functions  for  their  output  are 
provided.  These  will  be  discussed  in  the  next  chapter. 

3.3  Character  Maps 

Maps,  i.e.,  matrices,  whose  elements  are  characters 
(literals),  are  called  character  maps.  The  rules  stated 
for  value  maps  apply  also  to  character  maps,  except,  of 
course,  that  no  arithmetic  can  be  performed  on  the  latter. 
Information  can  also  be  extracted  from  them,  however,  by 
using  the  operator  "="  . 
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These  maps  are  widely  used  as  output  maps. 

Functions  for  character  map  output  are  included  in  the 
package.  Characters  that  may  be  used  include  all  valid 
APL  characters  including  some  overstrikes,  e.g.,  "$"  . 

Different  characters  are  used  to  denote  levels 
of  grey  in  output  maps.  Because  of  the  hardware  character¬ 
istics  of  the  terminal  and  the  very  low  printing  speed, 
overprinting  to  represent  shades  of  grey  is  not  desirable. 

As  can  be  seen  in  examples  using  characters  to  represent 
levels  of  grey,  characters  indeed  give  sufficient  resolu¬ 
tion.  There  can  be  any  number  of  grey  levels  in  the  same 
map  provided  that  there  are  enough  different  characters  to 
represent  them.  The  user  can  easily  change  the  characters 
in  a  character  map  and  get  the  output  represented  by 
different  characters,  until  the  visual  impression  is 
satisfactory  and  pleasing  to  the  eye. 

3. 4  Composite  Map 

Composite  maps  are  the  results  of  operations  and 
manipulations  performed  on  base  maps,  value  maps,  or 
character  maps.  The  ultimate  goal  of  this  package  is  to 
produce  composite  maps  by  combining  and  operating  upon 
given  base  maps  and  other  composite  maps. 

There  is  one  important  map  that  has  to  be  defined 
in  the  workspace  prior  to  any  operation.  It  is  used  to  accept 
the  results  of  operations,  and  contains  the  final  composite 
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map  to  be  output.  This  map,  the  working-map,  is  a  matrix 
named  MAP.  This  matrix  can  be  defined  using  the  APL  assign¬ 
ment  function  and  restructure  functions,  or  by  using  the 
DEFBWM  or  DEFCWM  functions  of  the  package. 

MAP  is  thus  a  reserved  word,  and  cannot  be  used 
for  any  other  purpose.  The  user  can  access  this  working- 
map,  change  it  completely,  change  particular  elements, 
i.e.,  points,  or  erase  it.  It  is  the  user's  responsibility, 
when  working  with  the  functions  of  this  package,  to  be 
careful  that  he  not  destroy  the  results  of  prior  operations 
which  are  in  the  working-map,  and  to  clear  it  before  his 
next  composition. 

If  the  working-map  is  undefined,  an  error  message 
will  usually  result.  This  will  also  occur  if  a  character 
or  binary  working  map  is  improperly  used.  The  functions 
for  output  of  composite  maps  are  incorporated  in  the  package, 
and  will  be  discussed  in  the  next  chapter. 

3. 5  APL  Functions 

Most  of  the  APL  functions  can  be  used  with  base 
maps  in  expressions  in  order  to  obtain  new  composite  maps. 
Arithmetic  dyadic  functions  are: 

+  addition, 

subtraction, 
t  division, 

x  multiplication,  and 

*  exponentiation. 
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When  using  these  functions,  APL  rules  apply.  For  example, 
when  a  matrix  is  multiplied  by  another  matrix,  the  result 
is  the  multiplication  of  corresponding  elements.  When  a 
matrix  is  multiplied  by  a  scalar,  the  result  is  a  multipli¬ 
cation  of  each  of  the  matrix  elements  by  that  scalar. 

Relational  and  logical  dyadic  functions  available 

are : 

<  less  than, 

<  less  than  or  equal  to, 

=  equal  to, 

>  greater  than, 

>  greater  than  or  equal  to, 

1  not  equal  to, 

a  and , 

v  or , 

a  not  and  (nand) , 

v-  not  or  (nor)  . 

The  logical  function  ~  (not),  which  is  monadic,  is  also 
available. 

These  functions  applied  on  base  maps,  or  other 
binary  maps,  will  give  binary  maps  as  a  result.  A  complex 
expression  can  be  made  by  concatenating  several  functions. 
Relational  and  logical  functions  are  performed  on 
corresponding  elements  of  two  matrices,  yielding  1  for 
true  or  0  for  false.  Combining  these  functions  with 
other  functions  in  an  expression  is  also  permitted, 
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provided  that  relational  and  logical  functions  always 
act  on  binary  maps. 


Whenever  any  syntactic  or  semantic  error  is 
encountered,  the  appropriate  error  message  is  printed 
and  execution  is  suppressed. 

Other  dyadic  functions  are: 


assignment  function, 
T  maximum  of, 

L  minimum  of, 
e  membership. 


The  assignment  function 
variable  name  a  value, 
point  of  the  map  can  be 
(MAP  [  1 ;  2]^-l)  will  assi 
the  working-map  in  the 


is  frequently  used.  It 
Using  APL  indexing  capab 
accessed.  For  example, 
gn  a  value  of  1  to  the 
first  row,  second  column. 


assigns  to  the 
ilities,  each 
the  expression 
element  of 


Monadic  APL  functions  are: 
x  signum, 

<j)  rotate  along  vertical  axis, 

0  rotate  along  horizontal  axis, 

®  natural  logarithm, 
absolute  value. 


There  are  also  other  APL  functions  that  may  be  used,  but 
they  either  restructure  the  matrix  or  require  its  prior 
restructuring.  The  result  of  applying  these  functions 
can  be  matrix,  vector,  or  scalar.  The  advanced  user 
will  find  out  how  to  use  them  by  consulting  APL  manuals. 
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APL  provides  very  efficient  tools  for  the  user  to 
create  a  variety  of  combination  maps  that  cannot  be  done 
by  hand.  Despite  the  relatively  slow  response,  the 
conversational  approach  makes  interaction  with  the  system 
as  convenient  as  possible.  The  functions  provided  by  APL 
PLUS  for  file  creation,  manipulation,  and  retrieval  are 
not  mentioned  here,  but  are  also  quite  easy  to  use. 

The  error  diagnostics  will  be  of  little  help 
to  the  user,  especially  when  the  programmed  functions  or 
the  nested  functions  return  an  error  message.  If  the  user 
does  not  know  APL  and  the  algorithm  used,  the  only  thing 
he  can  do  is  to  remove  the  suspension  of  the  erroneous 
function  and  try  again.  The  user  familiar  with  APL  and 
the  developed  functions  could  discover  within  the  function 
what  is  wrong  and  take  the  necessary  recovery  steps. 

The  same  applies  to  all  objects  in  the  workspace. 
The  user  can  examine  and  trace  out  the  necessary  information 
using  other  APL  functions  provided  for  that  purpose. 

3 . 6  Developed  Functions 

This  package  is  a  set  of  functions  programmed  in 
APL.  The  functions  can  be  used  independently  of  each  other, 
or  combined,  i.e.,  the  result  of  one  function  can  be  input 


to  another. 


-41- 


In  the  package  there  are  input  and  conversion 
functions,  area  filling  functions,  working-map  definition 
and  clearing  functions,  map  manipulation  functions,  contouring 
functions,  a  file-read  function,  map  generating  functions, 
and  output  functions.  Most  of  the  functions  are  syntactically 
niladic,  i.e.,  they  require  no  input.  Some  of  them,  however, 
need  parameters  in  order  to  do  calculations  or  make  decisions. 
Passing  parameters  to  these  functions  is  done  in  conversational 
mode,  rather  than  when  calling  the  function.  The  functions 
get  their  parameters  by  "asking”  questions  that  have  to  be 
answered.  If  a  wrong  parameter  is  input,  the  function  will 
either  ask  for  it  again  or  return  an  error  message.  The 
next  chapter  will  describe  each  of  the  functions  in  detail. 

All  functions  assume  that  there  is  defined  in 
the  workspace  a  working-map  named  MAP.  MAP  is  thus  a 
reserved  word  in  the  package.  The  dimension  of  the 
working-map  must  be  the  same  as  the  dimension  of  the  base 
maps . 

The  functions  can  handle  maps  of  any  dimension 
less  than  100  x  100  .  Actually,  this  limitation  is  imposed 
only  by  the  output  functions.  All  functions  assume  that  the 
directions  of  the  axes  are  as  shown  in  Figure  19.  Only 
positive  values  for  incides  of  the  matrices  are  allowed. 

If  coordinates  are  collected  on  the  basis  of  the  normal 
coordinate  system  (  y  upward) ,  they  must  be  transformed. 
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The  data  should  be  collected  on  the  basi 
rectangular  grid  used  on  IBM  2741  terminal.  This 
shown  on  Appendix  A.  In  this  case,  no  conversion 
The  use  of  this  particular  grid  was  determined  by 
hardware  because  the  print  positions  are  fixed. 


sofa 
grid  is 
is  needed, 
the 
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4.  FUNCTIONS  DEVELOPED  AND  ALGORITHMS  USED 

In  this  chapter  both  the  algorithms  and  the 
designed  functions  will  be  described.  The  algorithms 
are  APL-oriented.  Listings  of  the  APL  functions  are  given 
in  Appendix  D. 

4 . 1  Input  and  Conversion  Functions 

These  functions  are  designed  to  facilitate  the 
input  of  base  maps  and  their  conversion  from  a  square  grid 
to  rectangular  grid  points. 

4.1.1  MCONV  Function 

The  MCONV  (map  conversion)  function  is  used  for 
base  map  input  and  for  the  conversion  of  square  grid 
coordinates  to  rectangular  coordinates  when  so  specified. 
This  conversion  is  necessary  if  the  shape  and  area  of  the 
features  are  to  be  preserved.  Otherwise,  the  shape  will 
distort,  as  can  be  seen  in  Figure  9.  The  distorted  circle 
has  radius  11.  However,  since  11  units  in  the  y  direction 
are  longer  than  11  units  in  the  x  direction,  the  picture 
is  distorted.  On  the  same  figure,  the  undistorted  circle 
(generated  using  the  function  CIRCLE)  is  shown.  The 
distortion  has  been  removed  by  the  following  technique: 

The  ratio  between  x  and  y  units  is  7:11  .  Thus,  if 
all  x  components  are  multiplied  by  11/7  ,  the  shape  and 
area  of  the  figure  will  be  preserved.  As  the  multiplication 
gives  non-integral  values,  the  nearest  integer  is  chosen  by 
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simply  rounding,  i.e.,  0.5  is  added  and  the  smallest  integer 
in  that  value  is  taken.  X  =  smallest  integer  (xxll/7+0.5)  . 

The  whole  picture  is  then  displaced  so  that  the 
coordinate  with  a  minimal  value  for  x  is  the  same  as 
nonconverted  minimal  x  value.  Thus,  in  general,  the 
boundary  points  are  not  correct,  but  shape  and  area  are 
preserved  with  satisfactory  accuracy. 

Figure  10.  shows  the  square  defined  by  x 
coordinates  (10  15  15  10)  and  corresponding  y  coordinates 
(5  5  10  10)  both  before  and  after  conversion.  The  shape  of 
the  square  is  normal  after  conversion.  This  is  of  importance 
only  to  the  human  eye,  since  calculations  and  overlappings 
can  be  performed  with  distorted  maps,  with  even  greater 
accuracy  for  boundary  points . 

The  function  is  niladic  syntactically,  but  is 
designed  for  conversational  use,  and  thus  requires  that 
parameters  be  passed  to  it  in  answer  to  questions.  X  and 
Y  component  vectors  are  asked  for  first.  If  the  number  of 
elements  in  X  and  Y  does  not  correspond,  the  question 
is  repeated.  Next,  the  function  requires  information  about 
the  grid  to  which  the  coordinates  relate.  They  are  converted 
to  a  rectangular  grid  if  they  relate  to  a  square  grid.  When 
all  this  information  is  correctly  entered,  the  corresponding 
points  in  MAP  are  set  to  1. 
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Next  the  user  must  indicate  whether  the  points 
represent  an  open  or  closed  curve.  In  both  cases,  a  contour 
filling  routine  MCFILL  is  called  to  make  a  curve  of  the 
points  in  MAP. 

Optionally,  after  contour  filling,  map  output 
can  be  obtained,  and  X  and  Y  printed  out.  The  user 
can  see  which  points  define  the  contour.  This  is  helpful 
when  conversion  is  applied  or  when  there  has  been  a 
typing  error. 

In  case  of  an  input  error,  only  incorrect  X 
or  Y  elements  need  be  changed,  not  the  entire  vectors. 
MCONV  can  then  be  applied  again  to  the  new  vectors. 

If  the  points  entered  exceed  the  grid  boundaries 
(  X>70  or  Y>50  ),  they  will  eventually  be  ignored. 

However,  they  are  taken  into  account  when  MCFILL  is  applied. 
This  often  happens  for  points  collected  on  a  square  grid 
basis. 

If  the  boundary  of  the  grid  coincides  with  the 
boundary  of  the  area,  the  points  defining  the  boundary  line 
should  not  be  omitted  when  entering  X  or  Y  .  Only  two 
points  defining  that  straight  line  (the  first  and  last 
points)  need  be  entered. 

Again,  when  contour  filling  is  finished,  an  output 
map  could  be  obtained,  together  with  X  and  Y  . 
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4.1.2  MCFILL  Function 

This  function  is  designated  to  make  a  curve,  open 
or  closed,  from  the  set  of  points  given  by  vectors  X  and 
Y  .  This  function  is  automatically  called  by  the  MCONV 
function  whenever  the  user  enters  the  new  set  of  points 
representing  the  curve. 

It  can  be  used  independently,  but  X  and  Y 
must  be  defined  in  the  workspace  beforehand.  MCFILL  is 
also  a  niladic  function  that  requires  parameters  in  the 
conversational  mode.  Either  a  C  or  an  0  must  be 
entered,  depending  upon  whether  the  curve  is  closed  or 
open.  If  the  curve  is  specified  closed,  the  last  and  first 
pairs  of  coordinates  are  assumed  to  define  the  last  segment 
of  the  curve.  Otherwise,  the  curve  is  open,  and  this  segment 
is  not  included. 

The  function  assumes  further  that  the  curve 
segments  are  defined  by  points  in  the  order  in  which  they 
appear  in  vectors  X  and  Y  .  Thus,  subsequent  pairs 
X^,Y^  and  define  subsequent  segments  of  the 

curve  and  X2  2  define  the  first  segment; 

X2,Y^  and  X^jY^  the  second;  et  cetera).  Segments  are 
approximated  by  straight  lines. 

In  Figure  11.  the  starting  point  is  marked  1  . 
Subsequent  points  defining  straight  line  segments  are  shown 


as  small  circles. 
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The  algorithm  applied  in  this  function  assumes 
that  the  encoding  of  points  defining  the  enclosing  curve 
has  been  done  in  a  clockwise  direction.  In  Figure  11.  the 
arrows  indicate  the  direction  of  encoding  the  points.  The 
starting  point  is  irrelevant. 

If  there  are  two  points  at  opposite  sides  of 
the  ideal  contour  line,  the  algorithm  takes,  as  the  contour¬ 
line  point,  the  point  that  lies  inside  the  area. 

The  contour  line  filling  algorithm  is  outlined 
in  the  following. 

If  the  curve  is  designated  closed,  the  first 
component  of  each  vector  is  concatenated  to  its  end,  so  that 
Xn,Yn  and  1  (ecIual  to  X^,Y^  )  define  the  last 

segment  of  the  curve.  N  is  the  number  of  given  points  that 
define  the  curve. 

If  the  curve  is  designated  open,  X^,Y-^  are 
concatenated,  so  that  X^,Y^  and  X^,Y^  define  the  last 
segment,  which  obviously  has  length  0  .  In  other  words, 
there  is  no  segment  at  all. 

The  contour  line  filling  procedure  starts  with 
the  first  pair  of  points,  i.e.,  with  the  points  X(I),Y(I) 
and  X (1+1) , Y (1+ 1)  ,  where  1=1  .  If  they  are  neighbouring 

points  (each  point  has  eight  neighbours) ,  no  contouring  is 
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needed,  and  the  next  pair  is  considered  by  incrementing 
I  by  1  ,  and  starting  the  same  procedure.  Next, 

DX=X (1  +  1) -X  (I)  ,  and 

DY=Y (i  +  1) - Y ( I ) 

are  calculated,  and  corresponding  signs  are  stored  in  variables 
SX  and  SY  .  Depending  on  the  direction  of  the  line  segment, 
these  values  will  be  either  -1,  0,  or  +1  . 

Next,  vectors  XI  and  Y1  are  generated.  XI 
contains  all  possible  X  coordinates  that  contour  line 
points  could  have,  and  Y1  ,  all  possible  Y  coordinates. 

For  example,  for  the  line  segment  shown  in  Figure  12., 

XI  =(45)  ,  because,  looking  from  the  x  axis  between  the 

points  defining  the  segment  (x=3,y=7  and  x=6,y=3),  only 
two  points  can  be  included  on  the  segment  line,  viz.,  those 
with  x  coordinates  4  and  5  .  Y1  =  (6  5  4)  ,  because, 
looking  from  the  y  axis,  these  three  points  can  be  on  the 
line  segment.  Naturally,  3  points  should  be  contained  in 
the  ’'full"  line. 

If  | DX | > | DY |  ,  then,  for  each  point  of  XI  ,  the 

corresponding  Y  coordinate  is  calculated  and  stored  in 
the  vector  YP  .  The  following  formula  is  used: 

YP=Y(I)+ (Xl-X(I) ) xDY/DX 

If  | DY | > | DX |  ,  then,  for  the  points  of  Y1 

corresponding  X  coordinates  are  calculated  and  stored  in 
the  vector  XP  : 


XP=X(I)+tYl-Y(I) )  xDX/DY 
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Degenerate  cases  with  DX=0  or  DY=0  are 
handled  trivially. 

Thus,  in  the  example,  since  DX=3  and  DY=-4  , 
we  use  the  second  formula  to  compute  XP=(3.75  4.5  5.25)  . 

Because  XI,  Yl,  XP,  and  YP  are  all  vectors, 
their  components  are  all  given  their  values  "at  one  shot" 
in  APL.  However,  the  XP  and  YP  values  are  not  integral. 
Therefore,  prior  to  setting  the  corresponding  points  in 
MAP  to  1  ,  the  non-integral  coordinates  must  be  converted 
to  integers. 

As  can  be  seen  in  Figure  13.,  for  Yl(6  5  4)  , 

the  corresponding  XP  is  (3.75  4.5  5.25)  .  If  simple 
rounding  were  applied,  XP  would  result  in  (455)  . 

Simply  taking  minimal  or  maximal  integers  cannot  be  applied 
in  general.  The  point  with  X  value  =  5  and  Y  value  =  4 
lies  outside  the  area,  assuming  that  the  straight  line  is  an 
ideal  approximation  of  the  curve.  Taking  this  point  for  a 
contour  line  point  would  be  incorrect.  Therefore,  the 
following  algorithm  was  applied  to  include  as  contour  line 
points  only  those  that  lie  inside  the  area.  This  algorithm 
does  not  apply  to  open  curves.  In  the  case  of  an  open  curve, 
simple  rounding  is  quite  acceptable. 

For  each  coordinate  in  XP,  XPMIN  and  XPMAX  are 
calculated,  using  the  APL  functions  for  the  floor  (f)  of 
XP  and  the  ceiling  (L)  of  XP  .  Thus,  XPMIN  is  the  smallest 
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integer  in  XP  ,  and  XPMAX  the  greatest.  As  can  be  seen  in 
Figure  14.,  the  points  X(I),Y(I);  X (1+1) , Y (I +1)  ;  XPMIN,Y1; 
and  the  points  X(I),Y(I);  X  (1  + 1)  ,  Y  (1  +  1)  ;  XPMAX, Yl;  define 

two  triangles.  The  areas  of  the  triangles  can  be  calculated 
by  using  the  formula 


A  h[(X1-X2) (Y1+Y2J+(X2-X3) (y2+y3)+(X3-X1) (y3+Y1)] 

where  X^  z  X(I),Y(I)  ;  X2,Y2  -z  X(I*1),Y(W)  ;  and 
X3,Y3  =  XPMIN.Yl  for  the  first  triangle  or  E  XPMAX, Yl 

for  the  second.  The  areas  will  always  have  different  signs. 


In  the  normal  coordinate  system  (y  increasing 
upwards) ,  the  area  of  the  triangle  will  be  positive  if  the 
points  PI ,P2 ,P3  are  taken  in  counterclockwise  direction, 
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Figure  14. 


and  negative  if  taken  clockwise.  For  our  coordinate  system 
(y  increasing  downwards),  this  is  reversed,  so  that  the 
triangle  on  the  right  of  the  directed  line  segment  will 
have  a  positive  area,  and  the  triangle  on  the  left  negative. 
Because  of  the  encoding  rule,  all  the  points  on  the  right 
side  of  the  directed  line  segment  are  inside  the  area. 

Therefore,  the  point  giving  the  triangle  a  positive 
area  is  taken  as  the  contour  line  point,  despite  the  fact 
that  the  other  point  might  be  nearer.  For  any  point  which 
lies  on  the  segment,  the  area  of  both  triangles  is  0  , 
and  this  point  is  of  course  chosen. 


-54- 


This  procedure  is  applied  to 
segments.  The  corresponding  points  in 
set  to  1  .  If  the  points  lie  outside 
simply  ignored. 


all  points  of  all 
the  working-map  are 
the  grid,  they  are 


Finally,  the  (N+l)st  components  to  X  and 
are  dropped,  and  they  are  returned  to  their  original 
Figure  11.  shows  a  fairly  complex  curve  that  was  use 
test  the  function.  Input  points  are  marked  o  ,  and 
points  included  by  MCFILL  by  *  .  Figure  15.  shows 
same  area  in  a  usual  output  shape. 
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If  the  encoding  is  done  in  the  opposite  direction, 
viz. ,  counterclockwise,  the  function  performs  contouring  by 
taking  for  contour  points  only  those  points  that  lie  outside 
the  contour  area  or  on  the  segment. 

4 . 2  Area  Filling  Functions 

There  are  two  functions  to  fill  in  an  area  defined 
by  contour  lines:  MARSCAN  and  MARFILL.  MARSCAN  can  be 
successfully  applied  to  all  convex  areas  and  many  concave 
areas.  It  is  a  relatively  fast  function  compared  to  MARFILL, 
which  is  quite  general  but  very  time  consuming.  MARFILL 
examines  the  area  point  by  point,  and  tests  each  for  the 
many  various  cases  that  could  occur.  The  user  is  advised 
to  avoid  the  use  of  MARFILL  whenever  possible,  either  by 
using  MARSCAN  directly  or  by  dividing  the  area  into  a  set 
of  convex  areas  that  can  be  handled  by  MARSCAN.  These 
convex  areas  can  be  filled  by  MARSCAN  and  then  assembled 
together. 

4.2.1  MARSCAN  Function 

This  function  assumes  the  contour  of  the  area  to 
be  in  the  working-map  MAP.  The  first  part  of  MARSCAN  finds 
the  smallest  rectangle  that  contains  the  entire  area  to  be 
filled.  This  is  done  by  examining  each  row  until  the  first 
appearance  of  a  contour  point  is  detected.  The  same 
procedure  is  used  for  columns.  For  example,  to  find  MINX 
(minimal  X),  the  examination  is  started  on  the  first  column, 


-56- 


and  then  progresses  on  to  the  last  column.  To  find  MAXX 
(maximal  X) ,  this  procedure  is  reversed. 

When  the  smallest  rectangle  has  been  found,  the 
function  scans  each  row,  point  by  point,  until  a  point  of 
the  contour  line  is  detected.  This  is  illustrated  in 
Figure  16.  Scanning  is  done  in  the  four  directions  indicated 
by  the  arrows. 


1  2 

1234567890123456789012345 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

1 

2 


|  -■*»  *  * 
k  * 

\~p-  * 

|e»  * 

I* 
k* 
i-*-  * 

1-^* 

(  — ►  * 

| - p*.  * 


-  ►  * 


★  - - - - - - 

★  *  — . 

*  * 

*  *  -<3 - 

*  ★  *  - 

*  *  — 
*  -a 


* 

*  * 

*  * 


H 


* 

* 

*  “«* 

★ 


1  2 

1234567890123456789012345 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

1 

2 


\ 


3  [ — - ►  *  * 

*  it  -O - 1 

3 

i  a*  * A  i 

A*  *  k 

4  | - ►  *  * 

*  *  - | 

4 

i  i  1 1  !  u*  *a,  |  i 

A*  *  A 

5  | - ►  * 

*  -*3 - | 

5 

i]  Jhi  1M'!  !]j 

I!MI 

¥*t 


j  ¥* 

*¥ 

it*  *▼ 

V* 

*  Y 

T*  *T  i 

?  * 

f*  *¥ 

T* 

*  | 

r  *  T 

* 

A* 

A* 

A* 

|  A* 

* 

* 

* 

A  * 


*  A* 

*  A  A* 


*’ii 

*? 
★ 

*  A 
A  \ 


Figure  16. 


Initally,  vectors  XP  and  YP  are  empty. 
Coordinates  of  all  points  of  the  rectangle  that  are  scanned, 
i.e.,  do  not  belong  to  the  area,  are  concatenated  to  the 
vectors.  After  that,  the  whole  rectangle  is  set  to  1  in 
MAP.  Next,  all  the  points  in  XP  and  YP  are  reset  to  0  . 
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Figure  17.  shows  that  the  test  area  has  been  filled 
in  my  MARSCAN.  Comparing  Figures  17.  and  15.,  you  will 
notice  the  appearance  of  certain  superfluous  points: 

8,11;  9,11;  10,11;  and  32,13.  These  points,  however,  do 
not  exist  when  area  filling  is  performed  using  MARFILL,  as 
can  be  seen  in  Figure  18. 
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Figure  18 
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4.2.2  MARFILL  Function 

This  function  is  designed  to  fill  in  the  area 
enclosed  by  a  complex  contour  line,  including  those  that 
cannot  be  filled  by  MARSCAN. 

MARFILL  is  a  niladic  function  which  assumes  that 
the  vectors  X  and  Y  are  stored  in  the  workspace. 

The  algorithm  takes  advantage  of  the  grid  matrix 
representation  and  the  theorem  which  states  that  a  point 
is  inside  an  area  if  there  is  an  odd  number  of  intersections 
between  the  straight  line  originating  in  that  point  in 
any  direction  and  the  curve  that  defines  the  area.  The 
point  is  outside  the  area  if  the  number  of  intersections 
is  even  or  zero.  Figure  19.  illustrates  this  theorem. 

The  simplest  way  to  implement  this  would  be  to  take  the 
vertical  line  (in  one  direction  only)  originating  in  that 
particular  point  and  to  count  the  number  of  intersections. 
The  same  applies  for  each  point  of  the  map.  This  is 
essentially  what  the  algorithm  does. 
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First,  the  minimal  rectangle  containing  the 
defined  area  is  taken.  Obviously,  points  outside  the 
rectangle  cannot  belong  to  the  area,  thus  a  lot  of  useless 
calculation  is  saved.  Each  point  of  the  rectangle  is  then 
examined. 

Vectors  X  and  Y  contain  successive  pairs  of 
points  defining  straight  line  segments.  Suppose  that  X 
and  Y  have  N  components  each,  i.e.,  N  points  define 
the  curve  enclosing  the  area.  Also  suppose  that  I,J  are 
the  coordinates  of  a  point  inside  the  rectangle,  and  that 
X(K),Y(K)  and  X (K+ 1) , Y (K+l)  define  the  Kth  straight 
line  segment. 

For  a  straight  line  originating  in  P  ,  a  line 
parallel  to  the  y-axis  is  taken  in  the  same  direction  as 
the  y-axis.  This  is  illustrated  in  Figure  20.  A  non¬ 
degenerate  case  will  be  considered. 

First,  we  note  that  the  line  can  only  intersect 
with  the  segment  if 

( (X(K)  <I)  a  (X  (K+l)  >1) )  v  ( (X  (K)  >1)  a  (X  (K+l)  <I)  )  . 

In  other  words,  if  the  point  (I,J)  lies  between  the  points 
defining  the  segment,  there  can  be  an  intersection.  But 
this  condition  is  not  a  sufficient  one,  as  we  can  see  with 
point  P2. 
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Next,  we  check  the  following  condition: 

(Y(K) >J) a(Y(K+1) >J)  . 

If  this  holds,  there  must  be  an  intersection. 

Finally,  we  compute  the  coordinate  of  intersection 
between  the  entire  vertical  ray  and  the  segment: 

YK  =  Y(K)  +  (Y(K+l)-Y(K))x(I-X(K))/(X(K+l)-X(K))  . 

If  YK  >  J  ,  the  intersection  counts.  This  has  to  be  done 
for  every  combination  of  points  in  the  rectangle  and  straight 


segment . 


However,  this  is  merely  the  simplest  case. 
There  are  many  degenerate  cases,  and  each  combination 
must  be  examined  for  all  of  them.  The  functions  can 
handle  most  of  them  correctly. 


The  points  on  the  straight  line  segment  are  not 
important,  since  they  have  already  been  filled  by  MCFILL. 

For  each  point,  and  for  each  intersection  that 
counts,  1  is  added  to  the  variable  SUM.  After  examining 
all  line  segments,  the  corresponding  point  in  working-map 
MAP  is  set  to  1  ,  if  the  value  of  the  variable  SUM  is 
odd,  and  left  as  zero  otherwise. 

The  degenerate  cases  can  be  seen  in  Figures  18. 
and  21.  Both  areas  are  correctly  filled  in  by  MARFILL. 
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Areas  containing  "holes”  cannot  be  filled  in  by 
MARFILL.  A  simple  way  to  solve  that  problem  is  suggested 
in  Chapter  5. 

Another  degenerate  case  that  cannot  be  solved 
by  MARFILL  is  an  area  with  a  "tail",  as  shown  in  Figure  22. 
However,  this  can  also  be  easily  solved  by  entering  and 
connecting  the  tail  to  the  area  after  the  same  area,  without 
the  tail,  has  been  filled.  This  particular  area,  however, 
can  be  correctly  filled  by  MARSCAN,  but  it  does  not  work 
in  general. 


Figure  22. 
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Combining  MARFILL  and  MARSCAN  is  also  a  good  way 
of  increasing  the  efficiency  and  generality  of  area-filling 
functions  in  use.  It  can  solve  nearly  any  area  filling 
problem. 

Real-life  problems,  however,  are  not  so 
complicated,  and  most  of  them  can  be  solved  by  using  only 
MARSCAN.  MARFILL  is  therefore  more  important  theoretically 
than  practically.  Because  of  the  many  degenerate  cases, 
and  the  necessity  for  point-by-point  examinations,  it  is 
very  time-consuming,  and  should  be  used  only  when  necessary. 

4 . 3  Working-map  Definition  and  Clearing  Functions 

Four  very  simple  functions  are  developed,  mainly 
to  give  the  user  information  about  the  status  of  the  working- 
map,  MAP,  and  to  enable  its  quick  redefinition. 

4.3.1  DEFCWM  and  DE.FBWM  Functions 

Both  functions  are  used  to  define  the  working-map, 
called  MAP.  They  are  monadic  functions.  The  input  parameter 
required  is  the  two  component  vector.  The  first  component 
defines  the  number  of  rows  in  MAP  and  the  second,  the  number 
of  columns,  i.e.,  the  y  and  x  dimensions  of  the  grid 
respectively.  A  message  is  returned  to  the  user  that  the 
working-map  of  desired  dimensions  has  been  defined.  This 
is  useful  if  base  maps  of  different  dimensions  are  used. 
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DEFCWM  defines  a  character  working-map,  and 
DEFBWM,  a  binary  one.  These  functions  are  always  used 
when  beginning  the  session. 

4.3.2  CEMPTY  and  MEMPTY  Functions 

These  two  functions  are  designed  in  the  same  manner 
as  DEFCWM  and  DEFBWM.  They  are  used  to  clear  the  working-map. 
This  is  done  by  simply  taking  the  dimensions  of  the  existing 
working-map  and  redefining  them.  A  message  is  returned  to 
inform  the  user  that  the  character  or  binary  working-map  is 
empty  and  is  thus  ready  to  accept  new  information. 

CEMPTY  redefines  the  working-map  as  a  character 
map,  and  MEMPTY  as  a  binary  map,  regardless  of  its  previous 
state . 

4. 4  MAP  Manipulation  Functions 

Two  functions,  ASSIGN  and  ASSIGN1,  transform 
the  composite  binary  map  resulting  from  the  evaluation  of 
an  APL  expression  to  a  character  map  suitable  for  output. 

In  other  words,  these  functions  set  the  points  in  the 
working-map  to  desired  characters.  The  points  set  are  the 
points  of  the  binary  map  that  result  from  evaluation  of 
the  expression. 
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4.4.1  ASSIGN  Function 

This  function  is  dyadic.  The  left  argument 
should  be  a  character  enclosed  in  quotes  (literal) .  The 
right  argument  can  be  any  APL  expression  that  results  in  a 
binary  matrix  with  dimensions  equal  to  those  of  MAP. 

Each  point  in  the  binary  matrix  is  examined.  If 
its  value  is  1  ,  the  corresponding  point  in  the  working-map 
is  examined.  If  the  point  in  the  working-map  is  blank, 
the  desired  character  is  stored  in  that  position.  If  the 
point  is  not  blank,  its  coordinates  are  concatenated  to 
vectors  XP  and  YP  . 

If  all  the  points  have  been  examined  and  XP 
and  YP  are  not  zero-component  vectors,  the  message  is 
printed  to  inform  the  user  that  the  designated  points  have 
already  been  assigned  to  some  other  area(s).  The  coordinates 
(vectors  XP  and  YP  )  are  printed  out  as  well. 

The  user  can  specify  that  the  latest  assign  take 
place  or  that  it  not.  This  means  that  a  previous  assignment 
at  those  points  is  lost  or  preserved.  Thus,  the  ASSIGN 
function  should  be  used  whenever  there  is  a  possibility  that 
an  unwanted  overlapping  will  occur. 

4.4.2  ASSIGN1  Function 

This  function  does  the  same  as  ASSIGN,  except  it 
does  not  take  care  of  the  status  of  points  in  the  working-map. 
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The  new  assignment  takes  place  whether  or  not 
the  designated  points  belong  to  another  area.  It  does 
not  inform  the  user  about  any  overlapping. 

It  is  desirable  to  have  this  function  because 
there  are  sometimes  foreground  areas  and  background  areas 
on  the  map. 

If  the  ASSIGN  or  ASSIGN1  function  is  used  with 
a  value  map,  the  desired  character  is  assigned  to  all  the 
points  in  the  working-map  corresponding  to  non-zero  points 
in  the  value  map.  The  functions  cannot  be  used  directly 
on  character  maps. 

4. 5  Contouring  Functions 

Sometimes,  it  is  desirable  to  represent  an  area 
only  by  its  contour  line  in  order  to  achieve  better 
resolution  of  the  picture.  Three  functions  are  used  to 
obtain  the  contour  line  of  an  area  in  a  computer  map: 
CONLINE ,  CONLINE 1 ,  CONBLANK. 

These  functions  also  assume  that  the  character 
working-map  is  defined  in  the  workspace.  The  desired 
character  is  assigned  to  all  points  in  the  working-map 
that  correspond  to  the  points  of  the  contour  line  of  an 
area  in  the  binary  map. 
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The  functions  are  dyadic.  The  left  argument  is 
a  character,  enclosed  in  quotes,  that  should  be  used  as  a 
contour  line  character.  An  exception  is  CONBLANK,  where 
the  contour  is  blank,  and  the  remaining  area  is  represented 
by  the  character  entered. 

The  right  argument  can  be  any  valid  APL 
expression  which  results  in  a  binary  or  value  map  with 
appropriate  dimensions. 

4.5.1  CONLINE  Function 

This  function  is  based  on  the  fact  that  all  the 
points  which  do  not  have  a  complete  four-point  neighbourhood 
lie  on  the  contour  line.  The  four-point  neighbourhood  of 
X(I,J)  is  (I-1,J),  (1  +  1, J),  (I,J-1),  (I,J+1)  .  There  are 
problems  at  boundary  points  of  the  working-map,  because  a 
full  neighbourhood  does  not  exist  for  these  points.  The 
non-existent  points,  of  course,  cannot  change  the  status 
of  a  boundary  point.  For  example,  if,  for  a  boundary  point, 
all  the  remaining  neighbours  exist,  i.e.,  the  value  of  these 
points  is  unequal  to  0  ,  that  particular  point  is  assumed 
not  to  be  on  the  contour  line.  Otherwise,  it  is  included. 

To  be  able  to  apply  the  algorithm  uniformly  to  all 
the  points,  including  the  boundaries,  the  following  was  done. 
The  matrix  is  expanded  in  all  four  directions.  The  first 
row  of  the  matrix  is  concatenated  to  the  matrix  in  such  a 
way  that  it  becomes  the  first  row,  the  previous  first  row 
becomes  the  second,  and  so  on.  The  matrix  is  thus  expanded 
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in  one  direction,  and  its  first  and  second  rows  are  the  same. 
The  same  procedure  is  applied  to  the  last  row,  and  then  to 
the  first  and  last  columns.  The  algorithm  starts  examining 
the  second  element  in  the  second  row,  and  finishes  with  the 
(M-l)th  element  in  the  (N-l)th  row,  where  N  is  the  number 
of  rows,  and  M  the  number  of  columns  of  the  new  matrix. 

Because  the  matrix  was  expanded  in  the  described 
way,  the  new  rows  and  columns  do  not  affect  the  decision  as 
to  whether  or  not  the  point  lies  on  the  contour  line.  They 
do  allow  the  algorithm  to  be  applied  uniformly  to  all  points 
of  the  old  map  without  treating  boundary  points  as  a  special 
case . 

The  function  also  works  on  value  maps,  but  it 
is  slower  for  the  following  reasons.  Each  point  of  the  map 
is  examined.  If  its  value  is  equal  to  zero,  the  next  point 
is  taken.  Otherwise,  the  four-point  neighbourhood  is 
examined.  In  the  case  of  value-maps,  if  all  the  values  of 
the  four  neighbours  are  equal  to  the  value  of  the  point 
being  examined,  the  corresponding  point  in  the  working-map 
is  set  to  blank.  Otherwise,  the  point  is  set  to  the 
character  used  for  contouring. 

Thus,  in  a  binary  map,  points  with  value  0  are 
automatically  skipped,  and  only  the  relevant  points  are 
examined.  In  value  maps,  however,  all  the  points  with 
values  unequal  to  zero  are  examined.  Contours  are  obtained 
for  all  the  values  contained  in  the  map,  and  not  only  for 
one  desired  value. 
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When  all  points  of  the  unexpanded  map  have  been 
examined,  the  first  and  last  rows  and  columns  are  dropped; 
thus  the  input  map  is  returned  unexpanded,  hence  unchanged. 

Two  simple  examples  are  shown  in  Figure  23. 

Each  was  constructed  by  superimposing  three  base  maps  to 
which  CONLINE  had  been  applied.  Dots  are  assigned  to  the 
points  on  the  map  that  do  not  belong  to  any  of  the  areas 
shown. 
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4.5.2  C0NLINE1  Function 

This  function  is  essentially  the  same  as  CONLINE. 
However,  it  is  not  applicable  to  value  maps.  Also,  boundary 
points  of  the  map  (matrix)  will  automatically  be  taken  as 
contour  points  if  they  belong  to  the  area  being  contoured, 
as  well,  of  course,  as  if  they  really  do  lie  on  the  contour. 
The  reason  for  this  feature  is  that  it  makes  it  possible 
to  output  unframed  maps  which  are  sometimes  more  graphic 
and  eye-pleasing  than  a  standard  framed  map.  Simple 
illustration  of  that  is  shown  in  Figure  24.,  which  shows, 
in  a  new  shape,  the  same  two  examples  as  Figure  23.  Other 
examples  can  be  found  in  4.10. 
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To  achieve  this  contouring,  the  algorithm  in 
CONLINE  was  changed  only  in  that,  instead  of  concatenating 
the  first  row  to  the  matrix  itself,  the  logical  negation  of 
the  first  row  was  concatenated.  The  same  was  done  for  the 
first  column  and  the  last  row  and  column. 

The  effect  is  that  boundary  points  belonging  to 
the  area  shown  on  that  map  can  never  have  all  four  neighbours 
equal  to  1  ,  since  the  neighbours  in  the  expanded  row  or 
column  are  always  their  logical  negations.  This  function 
cannot  be  applied  to  value  maps  because  a  logical  negation 
of  an  integral  value  does  not  exist  in  APL. 

4.5.3  CONBLANK  Function 

Everything  stated  for  CONLINE  also  applies  to 
CONBLANK.  The  only  difference  is  that,  in  the  case  when 
all  the  values  of  the  four  neighbours  are  equal  to  the  value 
of  the  point  being  examined,  the  corresponding  point  in 
the  character  working-map  is  set  to  the  character  desired, 
and  is  blank  otherwise.  A  point  with  zero  value  is 
automatically  skipped. 

This  actually  results  in  a  reversed  CONLINE 
function.  Blanks  are  used  as  contour  line  characters, 
and  the  inside  of  the  area  is  represented  by  the  desired 
character.  The  two  examples  from  Figure  23.  can  be  seen 
again  in  Figure  25.,  after  being  operated  on  by  CONBLANK. 
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Figure  25. 

4 . 6  Map  Generation  Functions 

Three  functions  are  included  in  this  group: 

CIRCLE,  RECT,  and  STAR.  The  functions  are  used  to  generate 
the  geometrical  figures  indicated.  These  functions  may  be 
used  in  expressions,  provided  that  the  workspace  is  big 
enough. 

4.6.1  CIRCLE  Function 

This  function  is  monadic,  and  returns  an  explicit 
result.  The  right  argument  is  a  vector  consisting  of  three 
elements.  The  first  is  the  x  coordinate  of  the  centre  of 
the  circle  to  be  generated;  the  second  is  the 


y  coordinate; 
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and  the  third  is  the  radius  in  units  along  the  x  axis. 

Distorition  is  compensated  for  by  a  multiplication 
of  coordinates  by  the  distortion  factor  Q=ll/7  .  The 
function  assumes  that  the  working-map  is  defined  in  the 
workspace.  Its  dimensions  are  taken,  and  the  new  map  is 
generated  with  the  circle  centred  at  the  desired  position. 

To  generate  the  circle,  the  efficient  APL  function 

called  "inner  product"  is  used.  What  it  does,  essentially, 

is  to  examine  all  possible  combinations  of  x  and  y 

2  2  2 

coordinates  which  satisfy  the  relation  x  +  y  <  r  and 
set  their  values  to  1  .  Thus,  it  actually  sets  to  1 
all  points  that  are  inside  or  on  the  desired  circle. 

The  function,  because  of  the  use  of  the  inner 
product,  is  efficient  in  terms  of  speed,  but  requires  a 
lot  of  storage.  For  a  50  x  70  map,  it  is  nearly  impossible 
to  use  it  in  a  compound  expression. 

The  next  chapter  will  suggest  how  this  problem 
can  be  circumvented.  A  circle  generated  using  CIRCLE  is 
shown  in  Figure  9. 

4.6.2  STAR  Function 

This  function  works  in  the  same  way  as  CIRCLE, 

except  a  different  relation  is  to  be  satisfied,  viz., 

2_  _2_  2_ 

x3  +  y  3  <  r3  .  This  relation  yields  an  area  inside  the 
closed  curve  called  a  hypocycloid.  The  area  looks  like 
a  star  with  four  sides.  In  this  case,  the  distortion  is 


not  compensated  for. 
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This  function  was  used  to  generate  a  base  map  called 
HIGHWAY,  particularly  the  part  where  two  "highways"  cross 
each  other.  The  limitations  of  CIRCLE  also  apply  to  STAR 
because  of  the  inner  product  incorporated  therein.  Figure 
26.  shows  a  figure  generated  using  STAR. 
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Figure  26. 

4.6.3  RECT  Function 

This  function  generates  a  rectangle.  It  is  a 
dyadic  function  that  returns  an  explicit  result.  The  left 
argument  is  a  two-component  vector  which  specifies  the 
position  and  length  of  the  rectangle  in  the  x  direction, 
i.e.,  the  first  component  specifies  where  the  rectangle 
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starts,  and  the  second,  where  it  ends.  The  y  dimensions 
are  specified  by  the  right  argument  of  the  function. 

The  inner  product  is  used  to  take  all  the  pairs 
of  coordinates  that  satisfy  the  relation 

(  (Y<Y2)a(Y>Y1) ) a ( (X<X2)a (X>X1) ) 

and  to  generate  the  rectangle.  XI,  X2,  Yl,  Y2  are  the 
input  parameters  mentioned  above.  The  limitations  stated 
for  STAR  and  CIRCLE  also  apply  here.  The  Figure  27.  shows 
two  different  rectangles  generated  using  RECT. 
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4 • ^  Output  Functions 

Output  functions  exist  for  base-maps,  value  maps, 
and  character  maps.  All  of  the  functions  print  out  the 
working-map  MAP.  Therefore,  the  map  to  be  output  must  be 
placed  in  MAP  if  it  is  not  already  there. 

4.7.1  BMPRINT,  BMPRINT1  Functions 

BMPRINT  function  is  used  to  output  base  maps  and 
any  other  binary  maps.  The  frame  is  also  printed,  along 
with  the  numbers  that  indicate  the  coordinates  to  facilitate 
easier  orientation. 

The  function  is  niladic,  but  it  requires  some 
parameters  using  the  conversational  mode  of  interaction. 
First,  the  map  is  checked  for  values  of  zero  and  one.  If 
other  values  are  found,  the  message  "USE  VMPRINT"  is  printed, 
and  the  function  is  terminated. 

Next,  the  user  is  requested  to  enter  the  title. 

This  can  be  blank  or  a  string  of  up  to  120  characters,  i.e., 
one  line.  If  the  length  of  the  title  is  smaller  than  the 
length  of  the  map  in  the  x  direction,  the  title  is 
automatically  positioned  in  the  middle.  Otherwise,  it 
is  printed  starting  from  the  first  left  print  position. 

The  title  can  contain  any  valid  APL  character  except  the 


quote . 
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Two  characters  representing  two  levels  of  grey 
are  requested  next.  They  must  not  be  enclosed  in  quotes. 

The  first  character  is  used  for  0  values,  the  second  for 
values  of  1  .  A  blank  may  be  used  for  either. 

BMPRINT1  function  does  the  same  as  BMPRINT,  but 
no  parameters  are  requested.  It  is  used  for  quick  output. 
BMPRINT1  is  called  by  other  functions  such  as  MCONV,  MARSCAN, 
etc.  When  the  output  is  used  only  for  comparison  purposes 
or  as  an  auxiliary  output,  the  function  is  also  used.  The 
characters  used  for  representation  are  .  for  0  values, 
and  *  for  1  values.  Dot  is  used  because  it  makes  finding 
the  coordinates  of  the  point  easier  than  if  blanks  were  used. 
An  example  of  BMPRINT1  is  shown  in  Figure  28. 
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4.7.2  VMPRINT ,  VMPRINT2,  VMPRINT 3 ,  VMPRINT 4  Functions 

These  functions  are  used  for  value  map  output. 

The  user  is  actually  concerned  only  with  the  function 
VMPRINT.  The  other  functions  are  used  only  for  particular 
options,  and  are  automatically  called  by  VMPRINT  whenever 
respective  options  are  specified. 

First,  the  minimum  and  maximum  values  found  in  the 
map  are  printed  in  order  to  inform  the  user  what  values  he 
can  possibly  specify  to  be  represented  on  the  output  map. 

The  user  has  further  to  take  his  choice  of  four 
options,  and  to  specify  its  number.  The  options  are: 

1)  bare  map  -  a  map  without  any  description, 

2)  framed  map  -  a  map  with  a  frame  around  it, 

3)  framed  coordinate  map  -  a  map  with  a  frame  and 

numbers  to  indicate  the  coordinates  of  the  points , 

4)  coordinate  bare  map  -  a  bare  map  with  numbers  to 

indicate  the  coordinates  of  the  points. 

All  four  options  are  used  in  examples  in  4.10. 

They  are  also  included  in  CHPRINT  function. 

Further,  the  number  of  value  levels  and  the 
characters  representing  each  value  level  are  requested. 

A  value  level  represents  a  certain  value  range,  for  example, 
all  points  with  values  from  10  to  100. 
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From  each  level,  the  lower  and  upper  bounds  are 
requested.  No  two  levels  may  overlap,  i.e.,  no  value  can 
be  specified  to  belong  to  two  or  more  levels.  If  a  level 
contains  only  one  value,  this  value  must  be  entered  as  both 
upper  and  lower  bound. 

If  an  error  in  these  specifications  is  detected, 
an  error  message  is  printed,  and  the  specifications  are 
requested  again. 

Next,  the  title  should  be  entered.  The  map  is 

printed  row  by  row.  Values  that  do  not  belong  to  any  of 

the  specifications  are  printed  out  as  blanks. 

Figure  29a.  shows  a  small  value  map  printed  as 
a  matrix,  i.e.,  all  the  values  of  the  points  of  the  map. 
Figure  29b.  is  the  VMPRINT  output  of  this  map.  Only  four 

levels  were  specified.  The  first  level,  shown  by  *  , 
represents  all  values  on  the  map  from  1  to  12  ;  the 
second  level,  shown  by  □  ,  all  values  from  18  to  23  ; 
the  third  value  level,  shown  by  o  ,  all  values  from  45 
to  70  ;  and  the  fourth  level,  shown  by  |  ,  all  values 

from  90  to  100  .  Figure  29c.  shows  the  same  map  with 
the  bare  map  options.  The  first  level,  shown  by  *  ,  is 

from  1  to  12  ,  and  the  second  is  from  45  to  100  . 

Figure  30.  shows  a  bigger  map  generated  in  the  same  way  as 
the  small  one  of  Figure  27.  Values  that  do  not  belong  to 
any  of  the  four  levels  are  blanked.  Instead  of  graphic 
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Figure  30. 

character,  letters  or  numbers  can  be  used  for  level  repres¬ 
entation,  as  well  as  with  other  output  functions. 

4.7.3  CHPRINT  Function 

This  is  the  most  frequently  used,  and  the  fastest, 
output  function.  The  four  options  mentioned  in  4.7.2  are 
included  in  this  function.  This  function  is  also  niladic. 
The  user  is  requested  to  specify  only  the  title  and  the 
desired  option,  because  characters  are  already  contained 
in  the  map  itself.  The  map  is  printed  just  as  it  is, 
without  any  further  changes  or  conversions.  All  examples 
in  4.10  are  obtained  using  CHPRINT. 
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4. 8  File  Access  Functions 

The  only  function,  MREAD,  is  developed  more  for 
demonstration  purposes  than  for  real  use,  although  it  is 
employable  if  necessary. 

It  is  left  to  the  user’s  discretion  to  design 
the  file  system  according  to  his  needs.  More  will  be  said 
of  this  in  the  next  chapter. 

4.8.1  MREAD  Function 

This  function  is  monadic.  The  right  argument 
is  the  name  of  the  map  to  be  read.  It  must  be  enclosed  in 
quotes.  The  first  component  of  the  file  INDEX  is  read  in 
the  workspace.  INDEX  contains  the  names  of  the  maps  in 
the  order  in  which  they  are  stored  in  the  file  named  MAPS. 
Each  component  of  the  file  INDEX  is  compared  to  the  name 
of  the  map  to  be  retrieved,  until  the  correct  name  is 
found.  The  component  number  is  taken,  and  the  corresponding 
component  of  the  file  MAPS  is  read  in. 

The  function  can  be  used  in  an  APL  expression, 
but  the  retrieved  map  is  lost  after  the  expression  has  been 
evaluated.  To  keep  it  permanently  in  the  workspace, 
simple  assignment  should  be  used.  For  example, 

MAPONE  MREAD  ’RIVER’  ,  reads  a  map  called  RIVER  and  stores 
it  in  a  map  called  MAPONE.  There  are  many  other  solutions 
to  retrieve  the  desired  map  from  the  file. 
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4 . 9  Miscellaneous  Functions 
4.9.1  HISTOG  Function 

This  function  plots  a  histogram  of  the  values  of 
the  elements  of  the  input  vector  against  its  indices.  The 
function  is  dyadic.  The  right  argument  is  the  vector  with 
the  values  to  be  plotted.  The  left  argument  is  a  factor 
used  for  normalization  of  the  values,  i.e.,  the  constant 
by  which  all  values  are  divided. 

Each  value  is  first  rounded  by  dividing  by  the 
normalization  factor,  adding  0.5  ,  and  taking  the  smallest 
integer  in  that  value. 

HISTOG  can  plot  a  maximal  value  of  100  .  If  the 
values  to  be  plotted  are  greater  than  100,  a  message  is 
returned  to  use  an  appropriate  normalization  factor  to 
reduce  the  maximal  value  to  under  100.  The  function  could 
be  useful  in  exploring  visually  the  values  along  a  particular 
row  or  column  in  a  value  map.  The  examples  in  Figure  31. 
show  the  use  of  HISTOG.  MAP  was  plotted  along  its  29th 
column  and  22nd  row.  In  this  example,  MAP  contained  the 
sum  of  all  the  base  maps,  except  for  HIGHWAY.  The  same 
thing  can  be  seen  in  a  different  way  on  the  map  in  4.10 
showing  all  overlapping  points.  As  on  the  map,  the 
histogram  values  1,  2,  or  3  show  that  1,  2,  or  3 
maps  were  added  at  that  particular  point.  Another  example 
of  HISTOG  is  shown  in  Figure  32. 


CO  CM  T-i 
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Figure  31. 
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4.10  Examples 

Finally,  some  real  examples  will  be  shown.  All 
were  produced  using  only  seven  base  maps,  RIVER,  LAKE, 

SWAMP,  FOREST,  OVERBURDEN,  SLOPE,  and  HIGHWAY,  which 
appear  in  Appendices  B  and  C.  In  the  first  six  cases, 
points  defining  the  contour  lines  of  these  areas  were  input, 
then  the  MCONV,  MCFILL,  and  MARSCAN  functions  were  applied. 
HIGHWAY  was  generated  using  the  function  STAR.  In  order  to 
facilitate  the  comparison  of  the  base  maps  with  the 
composite  maps,  six  of  them  were  reduced  and  all  shown 
on  one  page  as  Figure  33. 

In  order  to  familiarize  the  reader  with  practical 
use  of  the  package  and  to  demonstrate  its  simplicity,  the 
entire  composition  work  will  be  shown  for  some  of  the 
produced  maps. 

The  computer  printout  in  Figure  34.  shows  the 
complete  session  to  compose  and  output  the  map  shown  on 
page  98.  Numbers  and  letters  at  the  beginning  of  each  line 
were  later  added  manually.  The  letter  0  indicates  that 
the  operator  at  the  terminal  (i.e.,  the  user)  typed  in  that 
particular  line;  the  letter  T  stands  for  the  line  printed 
by  the  terminal  as  an  answer  to  prior  input. 

Let  us  follow  each  particular  line  and  explain 
its  meaning.  In  line  1.,  the  user  signs  on  by  typing  his 
code  number  and  identification  key.  Lines  2., 


3.  and  4. 
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are  the  system’s  answer  containing  the  time  of  the  day, 
date,  and  users  name.  After  that  message  the  user  can 
start  his  own  work.  First,  at  line  5.,  the  library  called 
MAPS  is  loaded  into  the  workspace  (memory) .  Loading  is 
followed  by  printing  out  the  date  when  this  information  was 
saved.  Next,  on  line  7.,  the  functions  to  be  used  are 
copied  from  the  library  called  MAP FUN.  This  is  also 
followed  by  printing  the  saving  date  on  line  8.  Prior  to 
any  composition  work,  working-map  MAP  must  be  defined. 

DEFCWM  is  used  to  define  the  working-map  of  dimension 
50  x  70  .  This  is  an  example  of  a  monadic  function.  The 
system’s  acknowledgement  is  shown  on  line  10. 

The  real  composition  work  begins  on  line  11. 

The  dyadic  function  ASSIGN  is  used  to  represent  the  RIVER 
with  ”1".  On  lines  12.  and  13.  the  same  is  done  for  LAKE 
and  FOREST.  On  line  14.  an  APL  expression  is  used  as  a 
right  argument  of  the  ASSIGN1  function.  The  expression 
must  result  in  a  binary  matrix.  Points  in  the  working-map 
corresponding  to  elements  with  value  "1"  in  this  matrix 
are  set  to  the  desired  character.  On  line  15,  for  example, 
the  character  ”w"  is  used  for  all  points  in  the  SWAMP 
but  not  in  the  FOREST.  All  other  functions  that  accept 
APL  expressions  as  arguments  can  be  used  in  this  way. 

Line  16  shows  how  the  working-map  can  be  modified 
using  the  working-map  itself  in  the  expression.  Another 
way  to  do  the  same  thing  would  be 

’ . ’ ASSIGN1~ (RIVERvLAKE v FORE ST v SWAMP v SLOPE) 
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There  are  actually  many  different  ways  to  compose  the  same 
map.  Obviously  the  expression  on  line  16, is  much  shorter 
and  takes  less  time  for  the  system  to  calculate  it. 

The  user  unfamiliar  with  APL  notation  should 
bear  in  mind  that  APL  evaluates  expressions  from  right 
to  left. 

Actually  these  seven  lines  (10  to  16)  complete 
our  composition  work.  To  output  the  composed  map  CHPRINT 
is  used.  It  is  a  niladic  function,  hence  typing  its  name 
is  all  that  must  be  done.  Next,  on  line  18,  the  function 
asks  for  the  title.  It  is  entered  on  line  19. 

The  reader  should  not  be  confused  that  some  of 
the  lines  start  on  the  6th  printing  position.  Normally 
when  the  system  awaits  its  next  job  the  typing  ball  stops 
at  the  6th  print  position.  If  the  function  awaits  the 
answer,  however,  this  may  not  be  true. 

Next,  the  CHPRINT  function  requires  the  number 
of  the  desired  output  option.  It  also  informs  the  user 
that  the  new  page  should  be  positioned  (manually) ,  because 
immediately  after  the  RETURN  key  is  pressed  the  map  will  be 
output.  This  time  the  new  page  was  not  positioned,  so 
that  we  could  look  at  the  new  map  and  make  sure  that  it  is 
the  map  we  intended  to  compose. 

On  line  36,  the  ATTN  key  was  pressed  to  supress 
execution.  The  function  remains  suspended,  as  indicated 
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by  the  message  on  line  37.  To  remove  the  suspension, 

"->0"  was  typed  in.  To  end  the  session  the  user  signs-off 
by  typing  in  ")0FF". 

The  reader  can  see  that  the  composite  map  has 
been  created  in  a  very  simple  way. 

In  producing  the  next  16  maps  on  pages  99  to  114 
no  other  sophistication  was  used  than  that  explained  here. 
Instead  of  ASSIGN  and  ASSIGN1,  contouring  functions  were 
used  frequently.  To  compose  the  map  on  page  100,  for  example, 
the  following  functions  were  used  in  turn: 

C0NLINE1  RIVER 
C0NLINE1  LAKE 
’A’  CONLINE1  SLOPE 
'  oj'  CONLINE  SWAMP 
1 t '  ASSIGN  FORESTa- (SWAMP vSLOPE) 

Nonoverlapping  points  shown  on  page  113,  were  obtained 
using  only  two  instructions: 

MAP  l’e-SLO  PE  +  SWAMP +  RIVER+OVERBURDEN+ LAKE  +  FORE  ST 
f  0 ’ ASSIGN1  MAP1= 1 

All  maps  were  superimposed  and  the  resulting  map  stored 
in  matrix  MAPI. 

Overlapping  points,  shown  on  page  112,  were  then 
obtained  executing: 

CEMPTY 

’ 2 ’ ASSIGN1  MAP1=  2 
’ 3 ' ASSIGN1  MAP  1=  3 
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CEMPTY  was  used  to  clear  the  working-map.  The  number  2 
was  assigned  to  values  2  ,  i.e.,  the  points  where  two 
maps  overlap,  and  3  where  three  maps  overlap. 

In  composing  the  maps  on  pages  115-117  and  121-127 
the  CIRCLE  and  RECT  function  were  used.  These  functions 
utilize  a  lot  of  storage  and  cannot  be  used  in  the  expression 
because  WS  FULL  error  message  will  usually  result.  The 
following  example  shows  how  this  was  circumvented  (see  also 
Chapter  5) .  It  composes  the  map  on  page  122 


C0NLINE1  CHPRINT 


) CLEAR 

) COPY  MAPFUN  CIRCLE  RECT  DEFCWM 
DEFCWM  50  70 
CMAR^CIRCLE  54  24  32 
RMAP^5  40  RECT  3  10 
) COPY  MAPFUN  ASSIGN1 
) COPY  MAPS 
CONLINE1 
CONLINE 1 
CONLINE1 
A’CONLINEl 
o’ C0NLINE1 
* ’ ASSIGN1 
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RIVER 
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ASSIGN1 
ASSIGN 1 


RMAPa-(RIVERvSWAMPvHIGHWAYvOVERBURDEN) 
CMAPa~ ( LAKE v FOREST v SWAMP vSLO PE vH I GHWAY) 


CHPRINT 


The  workspace  is  first  cleared.  Only  CIRCLE,  RECT  and  DEFCW 
are  copied.  (The  answers  of  the  system  are  omitted  in  this 
description.)  Next  the  working-map  is  defined  and  the 
CIRCLE  function  used  to  generate  a  CMAP  matrix  that  contains 
all  the  points  within  the  circle  with  the  centre  in  X=54, 
Y=24  and  with  a  radius  of  32  X  units.  The  same  is  done  to 
generate  a  rectangle.  The  CMAP  and  RMAP  matrices  are  then 
further  used  in  logical  expression. 
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In  Chapter  5  some  other  suggestions  will  be 
given  concerning  this  and  other  frequently  occurring 
difficulties  and  problems. 

The  following  is  a  brief  description  and 
interpretation  of  the  sample  maps.  To  some,  a  title  was 
given;  often  it  was  deliberately  omitted.  Some  sequences 
illustrate  the  progressive  introduction  of  more  and  more 
information.  Others  illustrate  how  to  improve  the  visual 
quality  of  a  map.  One  such  technique  is  to  "separate” 
the  map  into  a  "foreground"  and  a  "background". 

The  first  example,  on  page  98  shows  a  composite 
map.  Space  containing  no  information  is  represented  by 
dots.  The  FOREST  is  shown  complete,  although  it  overlaps 
with  SLOPE  and  SWAMP.  This  can  be  interpreted  as  a  map 
showing  the  forest  and  other  surrounding  natural  environmental 
factors . 

The  next  example,  on  page  99  ,  shows  the  contour- 
blank  option  applied  to  the  previous  example. 

The  map  on  page  100  shows  the  contour  line  option  using 
the  function  CONLINE1.  The  output  option  is  called  a 
coordinate  bare  map.  SWAMP  and  SLOPE  are  shown  complete, 
i.e.,  the  map  shows  the  part  of  the  forest  that  is  neither 
on  the  slope  nor  in  the  swamp. 

,  The  next  map,  on  page  10 1,  shows  the  output 
option  called  a  framed  map.  Contours  and  full  areas 

can  be  combined  as  desired.  Different  characters  are  used 
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for  representation.  A  field  where  no  information  is  present 
is  shown  blank,  as  well  as  the  inside  of  some  areas. 

In  this  way,  three  levels  are  obtained  (contour, 
full  area,  blank);  these  make  a  good  visual  impression. 

The  full  areas  are  easily  noticed.  Thus,  this  way  of 
outputting  information  can  emphasize  some  items  in  the 
map . 

The  interpretation  of  that  map  would  be  the 
part  of  the  forest  that  is  not  on  the  slope,  and  the  swamp 
that  surrounds  the  forest. 

The  maps  on  pages  102  to  105  show  composite  maps 
of  four  base  maps.  Maps  are  output  and  produced  using 
different  output  options  and  composition  functions. 

The  next  map,  on  page  106,  shows  all  six  base 
maps  superimposed.  The  SWAMP,  RIVER,  LAKE,  and  SLOPE  are 
shown  in  full.  That  part  of  the  FOREST  which  is  not  on 
the  slope  or  in  the  swamp  is  shown  too.  OVERBURDEN  is 
seen  to  be  intersected  by  the  other  factors. 

To  improve  the  visual  impression,  the  CONTOURBLANK 
function  is  applied  to  the  FOREST  and  to  the  SLOPE  not  in 
the  FOREST.  As  can  be  seen  on  the  map  on  page  107,  two 
contour-blanked  factors  are  very  easily  noticed. 

The  map  on  page  108  shows  all  six  base  maps. 
Contoured  and  full  areas  are  combined.  It  shows  the  forest 
and  lake  as  foreground  information,  and  the  other  factors 
as  background. 
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The  next  three  maps  introduce  a  new  possibility. 

To  have  an  even  higher  contrast  between  the  foreground  and 
background  information,  the  contour  line  output  and  full 
areas  are  combined.  Additionally,  the  full  areas  are 
represented  by  an  overstrike  character.  The  complete 
forest  is  shown  in  the  first  of  these  examples.  The  part 
of  the  forest  which  is  neither  in  the  swamp  nor  on  the 
slope  is  shown  in  the  second,  and  the  complete  forest  with 
other  factors  on  the  third. 

The  map  on  page  112  shows  all  overlapping  points 
of  all  base  maps  except  HIGHWAY.  The  number  2  indicates 
that  there  is  an  overlap  of  two  areas,  and  the  number  3 
indicates  an  overlap  of  three  areas. 

Non- overlapping  points  are  shown  on  the  maps  on 
pages  113  and  114.  All  points  are  represented  by  the  same 
character  in  the  first  example.  In  the  second  example, 
points  belonging  to  different  areas  are  represented  by 
different  characters. 

Using  the  function  CIRCLE,  a  circle  is  generated. 

All  points  within  that  circle,  i.e.,  points  not  farther  than 
a  certain  length  from  a  certain  point,  and  not  in  the  swamp, 
on  the  lake,  or  on  the  slope,  are  chosen.  The  contour  blank 
function  is  then  applied.  Points  are  represented  by  an 
overstrike  character  to  emphasize  the  area  as  much  as  possible. 
This  map  is  on  page  115. 
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In  the  maps  on  the  next  two  pages,  two  other 
circles  are  taken  and  represented  by  overstrike  characters. 

On  contour  blank  output,  contouring  is  applied  to  some  of 
the  factors  (not  all),  but  not  to  the  circles.  As  can  be 
seen,  there  is  a  variety  of  possibilities. 

On  the  map  on  page  118,  the  HIGHWAY  is 
superimposed,  but  the  appearance  of  the  map  is  not  satisfactory 
because  the  picture  seems  a  little  vague.  On  the  next  map, 
only  one  character  representing  OVERBURDEN  is  changed.  The 
appearance  of  the  highway  is  then  much  better. 

The  same  map,  using  only  a  contoured  background, 
is  on  page  120. 

In  the  following  maps,  a  notion  of  possibilities 
for  planners  is  introduced. 

The  map  called  HOUSING  AREAS  on  page  121  is 
obtained  by  using  CIRCLE,  and  taking  all  the  points  in 
the  generated  circle  which  are  not  in  the  swamp  or  the 
forest  or  on  the  lake  or  the  slope  or  the  highway. 

"Industrial  area"  is  added  next  by  using  the 
function  RECT.  The  map  is  shown  on  page  122.  All  points 
in  the  rectangle,  but  not  in  the  river  or  swamp  or  on  the 
highway  or  overburden,  are  taken  and  represented  by  another 
overstrike  character.  Thus,  we  have  a  "housing  area"  and 
an  "industrial  area"  on  our  map. 
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On  page  123,  another  "project"  is  shown.  It  is 
completed  with  an  additional  "housing  area"  on  the  map  on 
page  124.  The  "industrial  area"  is  "planned",  however, 
regardless  of  the  other  factors,  so  that  it  takes  up  the 
full  rectangle.  Thus,  a  third  "small  project"  is  shown 
on  page  125. 

What  about  comparing  the  two  developed  projects? 
In  what  areas  do  the  "projects"  on  pages  122  and  123 
coincide  (overlap)?  They  are  very  easily  obtained  and 
displayed,  as  shown  on  page  126.  The  areas  of  "project" 
on  page  122  that  do  not  overlap  with  "project"  on  page  123 
are  shown  on  page  127. 
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5.  SOME  SUGGESTIONS  FOR  USE 

This  chapter  discusses  practical  aspects  useful 
for  persons  working  with  the  package. 

Input  of  areas  using  the  functions  of  4.1,  as 
already  mentioned,  cannot  cover  all  possible  cases.  For 
example,  areas  containing  a  "hole"  are  not  handled  properly. 
To  solve  this  problem  the  "hole"  should  be  defined  separately 
as  another  base  map  and  stored  in  the  workspace.  This  can 
be  done  by  using  the  assignment  operator 

HOLE  «-  MAP 

where  MAP  is  the  working-map  containing  the  area  that 
represents  the  "hole".  MAP  is  then  cleared  and  the  area 
that  contains  the  "hole"  defined  as  if  there  were  no  hole 
in  it.  The  base  map  MAPHOLE  that  contains  the  area  with 
the  "hole"  is  then  obtained  by  subtracting  HOLE  from  MAP 
and  storing  the  result  in  MAPHOLE,  e.g. 

MAPHOLE  MAP  -  HOLE 
or  MAPHOLE  +-  MAP  a~  HOLE  . 

Finally,  the  matrix  (map)  HOLE  is  erased. 

Using  that  method  of  restoring  intermediate 
results  and  combining  them,  many  of  the  problems  can  be 
solved.  Nearly  any  complex  area  can  be  obtained  by 
dividing  it  into  smaller  convex  areas  and  assembling  them 
together.  This  is  also  the  way  to  avoid  the  use  of  the 
MARFILL  function. 
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Many  of  the  intermediate  results  can  be  stored 
temporarily  in  the  user’s  library  and  retrieved  when  needed. 

The  library  can  restore  6  workspaces.  One  is  used  for 
storing  the  functions  of  the  package,  one  for  base  maps 
used  in  the  session,  and  the  other  4  are  free  for  the  user's 
intermediate  results.  For  example,  a  composite  map  frequently 
used  can  be  saved  in  this  way,  and  information  can  be  added 
to  it  or  deleted  to  produce  a  new  map  without  doing  the  entire 
composition  each  time. 

As  already  mentioned,  the  MARFILL  function  is 
very  time  consuming  because  of  point  by  point  examinations 
and  checking.  Therefore,  dividing  the  areas  into  convex 
polygons  is  a  much  cheaper  way.  The  faster  area  filling 
function  MARSCAN  can  be  applied  to  each  of  the  convex 
polygons  and  the  polygons  (areas)  assembled  together  to 
form  a  desired  map. 

If  the  points  are  collected  on  the  basis  of  the 
square  grid,  and  shape  is  to  be  preserved,  conversion  and 
"stretching"  the  area  will  take  place.  This  usually  results 
in  some  points  being  outside  the  grid.  The  whole  area  can 
be  shifted  back,  to  start  from  the  beginning  of  the  grid 
map.  This  can  be  done  by  giving  the  MCONV  function  the 
following  vector  as  input: 

X  «-  X  -  (L/X)  +  1 

This  could  actually  be  incorporated  as  an  option.  Unfortunately 
the  initial  coordinates  are  lost,  but  the  shape  and  integrity 
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of  the  area  are  preserved.  Initial  square  grid  coordinates 
are  lost  in  any  case,  because  otherwise  the  picture  is 
distorted.  Only  if  distortion  is  accepted  for  all  input 
base  maps,  the  coordinates  are  left  unchanged. 

The  use  of  the  package,  actually,  consists  of 
executing  the  functions.  Therefore,  most  of  the  errors 
occur  within  the  functions.  Whenever  such  an  error  occurs, 
the  user  should  type  in 

->  0 

to  remove  the  suspension  of  the  particular  function. 
Sometimes  this  should  be  typed  successively  two  or  even 
three  times  because  of  nested  functions. 

It  might  happen  especially  with  VMPRINT  function 
that  the  user  gets  caught  within  the  function  because  of 
improper  entry.  The  function  repeats  the  question  again 
and  again,  and  also  the  right  answer  does  not  help,  nor 
does  trying  to  escape  by  typing  -*  0  .  In  this  case,  the 
best  solution  is  to  type  in  the  nonsense  pattern  or  to 
continue  typing  in  wrong  entries  that  will  eventually 
result  in  an  APL  error  (for  example,  instead  of  a  vector 
the  user  can  enter  a  higher  dimensional  array,  which  will 
result  in  a  RANK  ERROR  message) .  When  the  APL  error 
message  is  printed,  typing  -*  0  will  bring  the  system  out 
of  the  function. 
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I£  the  user  is  experienced  and  familiar  with 
the  function,  much  information  can  be  obtained  by  displaying 
the  variables  of  the  erroneous  function.  This  should  be 
done,  of  course,  prior  to  removing  the  suspension,  because 
variables  local  to  the  function  disappear  when  the  suspen¬ 
sion  is  removed. 

The  functions  that  require  literal  input  should 
be  used  carefully.  If  the  user  gets  caught  by  such  a 
function,  typing  0  backspace  U  backspace  T  might 
help  him  to  get  out.  If  one  of  the  quotation  marks  is 
forgotten,  no  error  message  will  be  printed  because  APL 
accepts  everything  between  quotes  as  literal  data,  and 
as  long  as  the  closing  quotation  mark  is  missing,  no 
response  will  be  obtained  from  the  terminal. 

The  workspace  may  get  full  quickly  if  intermediate 
results  and  superfluous  functions  are  not  erased.  To  avoid 
this,  extensive  transfer  to  and  from  the  library  should  be 
used. 

The  storage  consuming  functions,  for  example, 
the  map  generating  functions  CIRCLE,  STAR,  RECT  should  be 
used  in  a  clear  workspace.  The  generated  map  should  be 
stored  as  an  intermediate  result  and  used  when  needed. 

These  functions  should  not  be  used  in  expressions  because 
the  WS  FULL  error  message  will  often  result. 
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Unfortunately ,  the  same  applies  for  value  maps 
if  real  numbers  are  used  as  elements.  Therefore,  the 
usage  of  value  maps  in  expressions  is  also  very  limited. 

Tricky  methods  and  combinations,  as  well  as  many  inter¬ 
mediate  steps,  must  be  used.  The  simplest  solution  and 
probably  the  best  one  at  this  point  would  be  simply  to 
request  more  storage  from  the  operator.  This  can  be  done 
using  either  the  terminal  itself  as  a  communication  device 
or  the  telephone  attached  to  it. 

Many  of  the  functions  accept  the  APL  expressions 
as  input.  For  example: 

ASSIGN1  RIVERv (SLOPEa-FOREST) a (MAP  =  '  ’) 

The  expression  can  be  of  any  complexity.  It  is  better  to 
split  the  complex  expression  into  a  few  less  complex  ones. 

The  part  of  the  expression  in  the  parenthesis  shows  how  the 
working-map  and  information  it  contains  might  be  used  to 
modify  the  working-map  itself.  In  the  same  way,  using 
relational  APL  functions,  information  can  be  extracted 
from  any  other  map  (value  map,  character  map). 

The  user  should  use  parentheses  whenever  appropriate. 
They  may  be  left  out  in  many  cases,  but  care  should  be  taken, 
especially  when  using  logical  functions,  because  of  the  APL 
syntactic  and  semantic  rules.  Remember  that  terms  of  an 
expression  are  associated  and  evaluated  from  right  to  left. 
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Some  of  the  functions  (MREAD, CIRCLE, etc. )  return 
an  explicit  result  when  used  explicitly.  It  is,  however, 
lost  when  the  function  execution  terminates.  To  retain 
the  result,  simple  assignment,  e.g.,  intermediate  storing, 
is  used.  For  example: 

CMAP  «-  CIRCLE  10  20  5 

will  store,  in  a  map  called  CMAP,  the  circle  with  centre  in 
X=10  Y=20  and  with  radius  of  5  units  (  X  axis  units) . 

CMAP  will  contain  0  values  for  points  outside  the  circle 
and  1  values  for  points  inside  the  circle.  CMAP  is 
used  in  further  manipulation  instead  of  CIRCLE  10  20  5  . 

The  user  should  further  take  care  to  define 
correctly  the  working-map  at  the  beginning  of  each  session, 
and  clear  it  when  doing  a  new  composition.  Some  functions 
require  a  character  working-map,  the  others,  binary  working-map. 
Inappropriate  definition  or  usage  will  result  in  an  APL  error 
message  DOMAIN  ERROR.  Using  map  definition  and  clearing 
functions  the  user  can  always  know  which  map  is  defined. 

When  output  functions  are  used  (VMPRINT,  BMPRINT, 
CHPRINT) ,  the  ATTN  key  can  be  pressed  to  stop  execution. 

The  function  is  suspended,  and  its  suspension  can  then  be 


removed. 
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6.  TIME,  COST  AND  STORAGE  CONSIDERATIONS 

It  is  always  of  interest  to  know  some  facts  about 
average  cost  associated  with  some  units  of  work.  Here  such 
an  estimation  will  be  given.  However,  since  it  is  impossible 
to  define  what  "average"  means,  an  estimate  will  be  given 
for  a  typical  session  in  which  one  map  is  composed  and 
printed. 

The  time  used  can  be  exactly  obtained  using 
APL  i  -  beam  functions,  that  derive  exact  CPU  and  connection¬ 
time.  Each  of  these  is  priced  separately.  For  example,  at 
the  University  of  Toronto  connection  time  costs  $3  per 
connection  hour  and  CPU  time  $3  per  minute.  These  prices 
will  be  used  in  the  examples. 

The  i  -  beam  function  i  19  (typed  in  as  an 
overstrike  character  consisting  of  t  and  i  characters) 
gives  the  total  time  the  user's  keyboard  has  been  unlocked, 
i  21  gives  the  CPU  time  since  sign-on.  The  time  is  given 
in  sixties  of  a  second.  To  have  the  time  printed  in  hours, 
minutes,  second  and  sixties  of  a  second 

60  60  60  60  ti 2 1 
and  60  60  60  60  ti19 

should  be  typed  in. 

The  workspace  that  takes  31,868  bytes  is  not  priced 
separately,  which  is  not  the  case  with  storage  for  files. 
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The  time  for  a  session  here  includes  time  spent 
on  all  necessary  manipulations,  e.g.,  loading  the  functions, 
defining  the  working-map,  applying  the  functions  to  produce 
the  composite  map,  etc.  In  the  example  here,  four  base 
maps  were  used  to  produce  a  composite  map.  The  C0NLINE1 
function  with  only  one  base  map  as  the  argument  was  used 
once,  ASSIGN1  twice,  and  ASSIGN1  once  on  a  logical  expression 
that  consisted  of  5  logical  operations  on  four  maps.  CHPRINT 
was  used  for  output.  The  same  applies  for  the  other  two 
examples.  Thus,  for  one  whole  session  producing  only  one 
map,  the  following  time  was  used  and  roughly  estimated: 

CPU  time  2  min.  27  sec.  $7.40 

Connection  time  13  min.  0.65 

Total  $8.05 

The  total  is  actually  computer  cost  for  producing  one  map  in 
one  session.  Time  for  only  printing  the  same  map  was: 

CPU  time  0.60  sec.  $0.03 

Connection  time  (m.s.)  0.33  (6  min.  30  sec.) 

Total  $0.36 

Time  spent  on  composing  the  map  (all  functions  and  base  maps 
already  in  workspace.) 

CPU  time  2  min.  24  sec.  $7.20 

Connection  time  5  min.  0.25 


Total 


$7.45 
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Naturally,  the  time  shown  is  very  much  dependent 
on  the  functions  used.  For  output  functions,  it  does  not 
differ  considerably,  but  for  manipulation,  contouring, 
map-generating  and  especially  area  filling  functions,  it 
differs  from  function  to  function.  It  also  depends  upon 
the  content  of  a  base  map  because  small  areas  are  handled 
faster  than  big  areas.  Finally,  it  depends  upon  the  computer 
system,  the  operating  system,  and  the  terminal  used  for 
output . 

Therefore,  these  numbers  can  be  considered  only 
as  an  indication  of  the  order  of  magnitude  that  can  be 
expected. 

In  the  examples  shown,  it  was  also  assumed  that 
the  workspace  was  big  enough  to  accommodate  the  base  maps, 
the  working-map,  and  the  overhead  resulting  from  expression 
evaluation.  If  this  were  not  true,  the  time,  e.g.,  the  cost 
increase,  would  be  unpredictable.  It  would  probably  be 
better  to  require  additional  workspace,  if  any  were 
available,  rather  than  try  to  circumvent  the  problems  using 
additional  and  tricky  manipulations  and  restoring  intermediate 
results.  However,  in  many  cases,  the  computer  system  is 
not  able  to  allocate  to  the  user  additional  workspace  upon 
his  demand.  Therefore,  in  many  points  throughout  this 
thesis,  it  has  been  suggested  how  to  solve  the  problem 
using  additional  manipulations,  e.g.,  time,  instead  of 
additional  storage  (workspace) . 
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7.  FUTURE  WORK  AND  IMPROVEMENTS 

Application  of  the  package  in  real-life  problems 
will  undoubtedly  suggest  some  changes.  However,  even  at 
this  point,  some  suggestions  for  possible  improvements 
and  for  further  work  can  be  made. 

Improvements  could  be  achieved  in  the  following: 
speed,  storage  utilization,  suitability  for  use,  error 
detection  and  error  messages.  Further  work  could  also  be 
done  in  developing  new  functions,  implementing  statistical 
and  mathematical  methods  used  in  computer  cartography, 
utilizing  the  force  of  APL  showing  the  new  possibilities 
and  suggesting  new  applications  in  various  fields,  linking 
the  package  to  other  I/O  devices  and  developing  general 
useful  APL  file  system  and  file-storing  concepts. 

Some  of  the  mentioned  improvements  and  further 
work  suggestions  will  now  be  discussed. 

Speed  is  probably  the  most  important  improvement 
needed  to  cut  the  total  computer  cost  per  map.  Here  are 
the  most  frequently  used  functions:  ASSIGN,  ASSIGN1, 
CONLINE ,  C0NLINE1 ,  CONBLANK.  They  should,  of  course,  be 
worked  on  first.  In  these  functions,  the  point  by  point 
approach  was  used  because  it  is  easier  and  solves  the 
problem  in  general.  If  this  approach  is  to  be  followed, 
the  first  improvement  that  could  save  time  significantly  is 
finding  the  minimal  rectangle  containing  the  area,  thus 
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avoiding  examination  of  all  points  in  the  map.  This  was 
used  in  area  filling  functions.  It  could  cut  the  CPU 
time  on  the  average  to  2/3  of  the  present  CPU  time. 

In  the  contouring  functions  the  improved  ASSIGN1 
function  could  be  used  to  set  all  the  points  of  an  area 
to  blank  (for  CONLINE  and  C0NLINE1  functions  only) .  After 
that  the  contour  line  of  the  area  should  be  followed,  and 
the  points  set  to  the  appropriate  character  or  blank.  This 
is,  however,  applicable  only  to  areas  enclosed  by  closed 
curves,  and  not  to  open  curves  or  discrete  points.  The 
present  functions  can  be  applied  to  discrete  points  or  open 
curves,  but,  because  they  do  not  have  a  contour- line ,  the 
result  will  be  blanks  (for  CONBLANK  function)  or  the  line 
or  points  itself  (for  CONLINE  and  C0NLINE1)  functions. 

Output  functions  are,  at  present,  satisfactory  in 
speed  and  suitability  with  one  exception.  VMPRINT  is 
impractical  in  that  the  user  can  get  caught  within  the 
function  and  cannot  terminate  the  execution  before  all 
answers  are  answered. 

Map  generation  functions  require  a  lot  of  storage, 
e.g.,  workspace,  but,  on  the  other  hand,  they  are  very  fast. 
Finding  the  smallest  rectangle  in  this  case  could  help  too, 
without  changing  the  speed  of  functions.  Functions  should 
generate  a  map  with  dimensions  of  the  rectangle,  and  then 
expand  it  to  the  normal  dimensions  of  the  base  maps. 
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Error  detection  and  diagnosis  could  be  incorporated 
in  some  functions.  According  to  the  errors  that  frequently 
occur,  some  algorithms  should  be  developed  and  implemented. 
Setting  some  counters  and  auxiliary  variables,  and  testing 
them  when  desired  and  needed,  would  be  a  better  approach 
than  implementing  clumsy  and  time-consuming  error  checking 
procedures,  error  messages,  and  error  recovery  procedures 
for  each  particular  case. 

Work  with  value  maps  is  a  serious  problem,  and 
could  impose  serious  limitations  to  the  use  of  the  package. 

It  is  not  caused  by  poor  design  of  the  package,  rather  by 
APL  and  the  size  of  the  workspace  allocated  to  the  user, 
on  the  one  hand,  and  by  the  dimensions  of  the  map,  on  the 
other.  Partitioning  value  maps  into  two  or  four  parts, 
storing  them  temporarily  in  the  user's  library, and  retrieving 
them  when  needed,  might  be  of  some  help.  The  speed  would 
be  considerably  decreased.  The  best  solution  is  to  double 
the  workspace  whenever  possible. 

Many  useful  functions  could  be  developed  in 
future  work: 

A  distance  function  could  answer  the  question: 
what  is  the  distance  between  two  points?  Area  and  length 
of  the  contour  determination  functions  might  be  of  use. 

Implementation  of  statistical  and  other 
mathematical  methods  used  in  computer  cartography  would 
add  new  application  possibilities. 
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A  function  to  divide  complex  areas  into  a  minimal 
set  of  convex  areas  is  needed  for  efficient  input  of  areas 
using  present  input  functions. 

A  function  to  inform  the  user  about  maximum  value, 
minimum  value  of  the  map,  number  of  areas,  size,  etc. 
could  be  developed. 

The  possibilities  of  APL  are  still  not  fully 
demonstrated  in  this  package.  APL  inner  and  outer  product 
functions  allow  the  APL  user  to  build  up,  in  a  very  simple 
way,  very  powerful  functions.  As  an  illustration,  the 
problem  involving  distances  between  the  cities  is  briefly 
mentioned  here.  The  distances  are  shown  in  a  form  of 
square  matrix.  They  are  not  necessarily  the  same  in  both 
directions  between  two  cities,  e.g.,  one  way  roads  may 
connect  two  cities  not  laid  out  parallel  to  each  other. 

If  the  matrix  is  called  D  ,  the  operation  DT.+D  will 
give  the  longest  two-leg  trips  from  any  city  to  any  other, 
and  DL.+D  will  give  the  shortest  two-leg  trip.  Instead 
of  T  or  L  functions,  any  other  appropriate  APL  function 
could  be  included,  thus  giving  another  meaning  to  the  inner 
product,  and  producing  another  matrix.  The  example  mentioned 
here  is  explained  in  (15.)  page  234.  The  outer  product 
function  is  similar. 

Obviously  this  is  a  simple  way  of  building  complex 
functions  that  can  be  applied  to  matrices,  e.g.  to  our  maps. 
Soon  there  will  be  available  APL  functions  allowing  access 
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to  a  printer  and  plotter  from  an  APL  terminal.  This  will 
open  the  door  to  many  other  possibilities  for  the  APL 
user . 

Data  bank  and  file-sharing  concepts  will  also  be 
important  in  the  future.  Further  work  should  definitely 
be  attempted.  The  development  of  a  powerful  and  versatile 
APL  package  for  computer  mapping,  or  even  an  APL-based 
language,  would  certainly  benefit  the  development  and 
acceptance  of  computer  mapping  on  a  wide  basis. 
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8.  CONCLUSION 

Finally,  here  is  a  word  of  conclusion  drawn  from 
experience  and  impressions  while  designing,  implementing, 
debugging,  and  testing  the  package,  and  \\rorking  on  examples. 

The  APL  language  and  its  comprehensive  and  powerful 
instructions  were  found  much  easier  to  use  than  expected, 
and  very  natural  for  expressing  derived  operations.  Its 
debugging  facilities  were  found  satisfactory. 

Its  relative  slowness  should  not  discourage  the 
user;  on  the  contrary,  it  allows  him  to  take  a  "breath" 
and  decide  about  the  next  step.  The  conversational  approach 
has  been  found  attractive  and  challenging,  although  not  as 
much  as  the  interactive  approach  using  a  CRT  might  be.  Unless 
the  cost  is  too  high,  there  should  be  no  reason  for  non- 
acceptance  of  APL  and  conversational  approach  in  computer 
cartography . 

Future  urban  geographical  and  similar  data  banks 
and  APL's  file-sharing  facilities  and  interactiveness  are 
other  arguments  for  its  acceptance. 

The  advantages  of  the  package  are  that  the  user 
in  a  fairly  simple  conversational  way  can  produce  maps, 
change  them  immediately,  and  get  the  hard  copy  output. 

He  can  do  it  when  he  pleases  without  any  intermedia  and 
without  any  other  persons  for  intermediate  operations 
(programmer,  keypuncher,  etc.). 
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The  disadvantages  are  the  need  for  a  relatively 
big  storage  (workspace)  whenever  full  capabilities  of  APL 
are  to  be  utilized,  and  the  relatively  small  mapping  area 
that  can  be  handled  and  printed. 

This  can,  however,  often  be  remedied  and  circumvented 
by  using  various  techniques  and  sophisticated  composition  of 
maps,  at  the  cost  of  increasing  the  time  needed  for  the 
particular  work. 

It  appeared  that  it  is  very  difficult  to  think 
in  terms  of  overlapping  areas.  Many  times  the  results  were 
different  than  imagined  and  expected,  although  the  complexity 
of  the  expressions  were  not  great.  This  indicates  too  that 
the  package  is  not  only  a  tool  for  reproducing  maps.  It 
is  rather  a  tool  for  composing,  creating  and  producing  new 
maps  from  existing  ones.  These  produced  maps  cannot  be 
easily  done  by  hand.  It  is  a  real  automated  mapping  tool 
although  still  with  limited  possibilities. 

Thus,  the  first  introduction  of  APL  and  a  remote 
terminal  into  the  field  of  computer  cartography  has  been 
done,  and,  I  hope,  successfully. 
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APPENDIX  A 
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RECTANGULAR  GRID  USED  ON  IBM  APL  TERMINAL  2741 
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APPENDIX  B 
Input  Base  Maps 
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APPENDIX  C 
Generated  Base  Map 
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APPENDIX  D 
Functions  printouts 
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V  MCONViAiXliYliIiXP 

Cl]  ' ENTER  VECTOR  OF  X  COMPONENTS' 

C  2]  X+U 

[3]  ' ENTER  VECTOR  OF  Y  COMPONENTS' 

[4]  Y+Q 

[5]  -*-Lixi(+/x=x)  =  (+/y=.n 

[6]  • INCORRECT  NUMBER  OF  X  OR  Y  COMPONENTS' 

C7]  -1 

[8]  LI i' ENTER  YES  IF  COORDINATES  RELATE  TO  SQUARE  GRID ' 

C  9 ]  -*\£2x  i  a/  (  ,y£,5,=3pi4-m!]) 

CIO]  X+U+0.5 
C  11  ]  -*L3 

C  1 2  ]  L2  :  XP+-\./ X 
C  1 3  ]  X-«-L(Xxll*7)  +  0.5 

C 1 4 ]  X+X-((l/X)-XP) 

C 15 ]  L3:Y+lY+0.5 
C 16 ]  1+1 

c  17]  m:xi+,xin 
c  i8  ]  yn-.yci] 

C 19 ]  +L5x\(Xl>~l*pMAP)v(Yl>lipMAP) 

C  20]  MAPZYliXll+l 
C  21  ]  L5  : 1+-I  +  1 
C  22]  -+L 4xi  (  J  s  p  X  ) 

C  23 ]  'ENTER  YES  IF  CONTOUR  FILLING  DESIRED ' 

C  24]  +L6xia/(  »yES’*3p/l+a) 

C  25 ]  MCFILL 

C  26]  L6:' ENTER  YES  IF  MAP  OUTPUT  DESIRED .POSITION  NEW  PAGE  PRESS  RET.' 
C  27  ]  -►Oxaa/C  •r£,S,*3pyl<-a) 

C  28 ]  BMPRINT 1 
C  29 ]  'X' 

C  30]  OX 
C  31 ]  'Y' 

C  32]  O Y 

V 
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V  MCFI LL  -.A  i  B  ■.  CUR  ■.  DX  DY ;  I-.J-.PX  -.PYiPl  -,P2 SX -,SI ;  XP-,  XI  lYP-.Yl 

[1]  i6:’ ENTER  C  IP  CONTOUR  IS  CLOSED  0  IF  IT  IS  OPEN  CURVE* 

[2]  CUR+ lpB 

[  3]  +m*\CUR=*  c* 

[4]  -►iS  *  \CUR-*  O' 

[5]  +L6.P  'ERROR  ■.  SPECIFY  AGAIN* 

[6]  i4  ( 1  tX) 

[7]  y-*-y,(ity) 

[8]  -*•£  7 

[9]  iS:y^y,("i+y) 

[10]  Y+Y,C  1  +  X) 

[11]  LT-.I+l 

[12]  £8:PX-U[X]=A'[x+i]+i)vU[x]=y[Xi-i]-i  )v(y[x]=y[x-n]) 

[13]  py*-(y[x]=y[x+i]+i )v(y[x]=r[x*i]-i)v(y[x]=y[x+i] ) 

[14]  +L\*\PX*PY 

[is]  sx+*oy-y[x+i]-/[j] 

Ci6]  sy»a£>r*y[x+i]-r[x] 

[17]  Xl*X[X]  +  S/*/l*i  I  (y[X  +  l]-.Y[X])-(  1*SX) 

[is]  yi*-x[X]+sy*B*-i  i  (y[x+i]-x[x])-(i *sy) 

[19]  •*£2*i(  |£>X)2(  \DY) 

[20]  XP*-(X[X]-(£X»£y)*X[X]  )+Xlx(DAr*£iX) 

[ 21 ]  yp-xi 

[22]  -*L  \2*  \CUR-  '  0  ' 

[23]  *L3 

[24]  i!2:XP*UP+0.S 

[25]  -»i  9 

[26]  I2iXF*~Xi 

[  27]  YP-> {YlI]-(DY*VX)*XlIi  )+Xl*(.DY*DX) 

[28]  -+L13  *  \  CUR=  *  0  ' 

[29]  -►i  1 4 

[30]  il3:yP*[XP+0. 5 

[31]  -»i9 

[32]  £3sPi-(i  +  2)x( (y[x]-y[x+i] )x(y[x]+y[x+i]))  +  ( (y[x+i]-[yp)*(y[x+i]+yp))  +  (( ([yp)-y[x])«(yp+y[x] ) ) 

[33]  p2-(i*2)*((y[x]-y[x+i])x(y[x]+y[x+i]))+((y[x+i]-ryp)x(y[x+i]+yp))+(((ryp)-y[x])*(yp+y[x])) 

[34]  XP*(.  (UP)x(Pl>0)  )♦(  (  ry?)x?2>0)t(  ([yP+0.5)*(Pl  =  0)vP2  =  0) 

[35]  *L 9 

[36]  ii4:Pi*(  i  +  2)x(-Dyxy[x]+y[x+n )+((y[x+i]-yp)«(y[x-n]+iyp)  )+((/p-y[X])x( (iyp)+y[X]) ) 

[37]  P2*(i  +  2)x(-Dyxy[j]+y[j+i] )  +  ( (x[x+i]-yp)x(y[x+i]+ry?) )  +  ( u?-x[x] ) x( < ryp) +y[x] ) ) 

[38]  yp*(  ([yp)xpi>o)  +  ( (ryp)xP2>o)+((iyp+o.5)x(pi  =  o)vP2  =  o) 

[39]  i 9:^-1 

[40]  iio:-»tii*i(yp[x]>“itpW/iP)v(yp[x]>itpWi4P) 

[41]  «4P[yp[/];yp[^]]-i 

[42]  ill  :+il0xi(  J+J  +  \  )<.pXP 

[43]  il:-»£8x»(X-<-X-t-l  )s(pX)-l 

[44] 

[45]  yx-'i+y 

7 
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7  MARSCAll  \MINX  \MINY  \MAXX  \MAXY  -,I  \J  -%L\MiXP  \YP 
Cl]  MINX+MINY+ 1 
[23  MAXX+~\\pMAP 

[3]  MAXY+lipHAP 

[4]  I+J  «-l 

[5]  LI :  -*\L2x  \  (  +/MAPIJ ;  ]*0 ) 

[6]  ■+Llx\(MINYW+J  +  l)<ZlipMAP 

C7]  -*0,pD  +-' THERE  IS  NO  AREA  IN  THIS  MAP ' 

[8]  L2:*7«-l+pM/lP 

[9]  L3:+L4xi (t/MlPC^; ]*o) 

CIO]  -►L3x(M/i;ry«-.W-l  )>0 
Cll]  L4:-*L5xt(  +/.A//1PC  ;J]*0) 

C  1 2]  -*-mx(MINX’*-I-*-I+l ) W pMAP 

C 13 ]  L5:I<-~lipMAP 
C 14]  LS;+L7xx(+/ 'MAPI ;I]*0) 

C  15  ]  -*L6  x  ( MAXX+-I+-I  -1 )  >  0 

C  16]  LI :  XP+-YP-*-\  0 
C 1 7 ]  M+MINX 
C 18 ]  L+MINY 

C  19]  LI  5  :-»-Ll  6x  iMlP[L  ;M]  =  1 
C20]  XP+XP,M 

[21]  YP+YP.L 

[22]  +L15xx(L+L+l)zMAXY 

[23]  L16-.L+MINY 

[24]  +L15xx(M+M+1)<;MAXX 

[25]  L+-MAXY 

[26]  M+MINX 

[27]  LI  7  :-*L18x\/L4P[L;W]  =  l 

[28]  XP+XP.M 

[29]  YP+YP ,L 

[30]  -*L1 7  x  \  ( L«-L  -I  K  MI  NY 

[31]  L1Q-.L+MAXY 

[32]  -*L17xi  (M+M+1)£MAXX 

[33]  L+MINY 

[34]  M+MINX 

[35]  LI 9  : -+L20 x i/L4P[L  ;M]  =  1 

[36]  XP-XP.M 

[373  yp+yp.l 

[38]  +L19X!  (M+M+l  )<.MAXX 

[39]  L20  :M*-MINX 

[40]  +L19xX(L*-L  +  l)zMAXY 

[41]  L+MINY 

[42]  M+MAXX 

[43]  L21 i+L22xxMAPLLiMl=l 

[44]  XP+XP.M 

[45]  yp^yp.L 

[46]  +L21xX(M+M-l)zMINX 

[47]  L22:M+MAXX 

[48]  +L21xx(L4-L  +  1)ZMAXY 

[49]  MAPlMINY+\(MAXY-MINY  +  l)\MINX+\{MAXX-MINX  +  \  )  >1 

[50]  I<*-1 

[51]  -+L  24x\  (  p  XP  )  =  ( p  \  0  ) 

[52]  L23:ALlP[yP[I]  ;*?[!]  ]<-0 

[53]  -*>L23xi  ( J-e-J  +  l  )<>pXP 

[54]  'ENTER  YES  IF  MAP  OUTPUT  DESIRED ' 

[55]  +0xx*/  '  YES'  3pP)) 

[56]  L24 : BMPRINT1 
V 
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7  MARFILL \YR\XR\YL\XL ;MINX ;MINY \MAXX \MAXY \SUM\I \J \K\A1\B1 \K\\lh%YK 
Cl]  N*-~ 1 

[2]  114:  -►LIU x  i  (  XR+1  t  (N*-N  +  l  )cj»AT )  =  XC  p^D 

[3]  Yfl^ltU^Y) 

[  4  ]  N*-  1 

C  5  ]  L15:+L15xi(  XL+~li  (  -N*-N  +  l)QX)  =Xll] 

[6]  YL+~lib  -//) 4>1) 

[7]  X«-(“lt;D,**(l  +  X) 

[8]  y4-("ity),y,(itjr) 

[9]  MINX+l/X 

[10]  MAXX<r\/X 

[11]  MINY+l/Y 

[12]  MAXY*-[/Y 

[13]  I+MINX 

[14]  19 iJ+MINY 

[15]  16 :K*-2 

[16]  SUM*-  0 

[17]  12:->llxi((y[X]<J)A(y[^+l]>J))v6  y[X]>I)A(y[^  +  l]<J)  ) 

[18]  -+L3*\XlKl*I 

[19]  -►112xiy[^]  =  y[X  +  l] 

[  20]  *►(13,11  )[i  +  (UU-i]<i)aU[*+i]>i)  )v(  (y[^-i]>j)A(y[/f+i]<j) )] 

[21]  112:  -►13  xxYlKl^J 

[22]  ->L16xltf<(py)-l 

[23]  Bl+XR 

[24]  Al+XlK- 1] 

[25]  Kl+K+1 

[26]  *>117 

[27]  116 : Al+XZK-l ] 

[28]  Kl*-K 

[29]  113:  -►11 3  x  i  AlKl+tfl  +1]=Y[K1] 

[30]  B1+XZK1] 

[31]  117:^14xl(  U1<I)a(B1>I)  )v(  U1>I)a(B1<I)  ) 

[32]  -►(15,13)[l  +  (  K*-Kl  -l)£(pY)-l] 

[33]  L3:+L2xi(K«-tf+l)s(p*)-l 

[34]  -►15 

[35]  11  :+L3xi  ( Y[ £]<*/)  a(Y[*  +  1]<.7) 

[36]  Y/^Y[/n+((YU+i]-jm)*uu+i]-Y[/n )  )x(i-*m) 

[37]  -►lUxt  YK'Z.J 

[38]  -+L3 

[39]  L*:SUM*-SUM+1 

[40]  -+L3 

[41]  15  :-^L8  x  i  0  =  2  I  SUM 

[42]  MAPlJ ;J]«-1 

[43]  18  :-^16xi  (*W  +  l  )*MAXY 

[44]  -►19 xi  ( J^-I  +  l  )$MAXX 

[45]  BMPRINT1 

[46]  X<-l+X 

[47]  X*-~l+X 

[48]  y«-i+y 

[49]  Y*-~l+Y 
V 


I 
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V  DEFBWM  RC 

[1]  MA P+RC p  0 

C2]  'BINARY  WORKING-MAP  *  ;/?£[!];'  BY  *  ;i?CC  2 3  ;  •  DEFINED ' 

V 


V  DEFCWM  RC 
Cl]  MAP+RCp '  ' 

[2]  • CHARACTER  WORKING-MAP  '  ;/?£[!];'  BY  '  \RCl  2];’ 


V 


DEFINED ' 
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7  CEMPTYiRiC 
Cl]  R+lipMAP 
[2]  C+~lipMAP 
C  33  MAP+(R,C)p'  ' 

[4]  * CHARACT .  WORKING-MAP  *;/?;'  •  s C ;  *  EMPTY' 

7 


7  MEMPTYiRiC 

[1] 

[2]  C+~l+pMAP 

[3]  MAP-*-(R  ,C  )  p  0 

[4]  • BINARY  WORKING-MAP  '  ; i? ;  *  Bl  '  ;C;'  EMPTY' 
7 
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V  CH  ASSIGN  EXPR ; CR\XP\YP\YES\I \J \K 
Cl]  CH+lpCH 

C  23  XP+YP+ iO 

C  33  J+1 

C4]  L5 i I +1 

[5]  L^:+Llx\EXPRZJ ;I]  =  0 

[6]  -+L2x\MAPlJ ;  J]*  *  * 

[7]  MAPlJil^+CH 

C  8  3  +L1 

C  93  L2iXP+XP,I 

[10]  YP+YP ,J 

[11]  Ll:-*-L4x\(j-^j+i  )z~l  +  pEXPR 

[12]  -*•£ 5*  \  (J  +-J + 1 )  £  1  +  p EXPR 

[13]  **-0x  i  (  pXP )  =  0 

[14]  'THIS  AREA  OVERLAPPS  WITH  AN  OTHER  AREA (S)  IN  POINTS 

[15]  *X- * i XP 

[16]  'Y-'i YP 

[17]  • ENTER  YES  IF  LATEST  ASSIGN  FOR  THESE  POINTS  TAKES  PLACE » 

[18]  y£S«-3pH 

[19]  -*-Oxia /'YES'*YES 

[20]  K+l 

[21]  L7  -.MAPlYPlKl  iXPlKll+CH 

[22]  -*-Ll*\  (K+K+l  )£pXP 

V 
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7  CH  ASSIGN 1  EXPR\I\J\CH 
[13  CH+1qCH 
[2]  J+ 1 
[33  L5:I«-1 

[4]  L  4  :  -*-L  1  x  i  EXPR  [  J ;  J  3  =  0 
[53  MAPUill+CH 
[63  Ll:**-Z/4xi(jVJ  +  l)5“l  ip  EXPR 
[73  **£5xl  (J+J  +  l)zlipEXPR 


7 
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V  CH  C0NLINE1  E \D  ;i?  iE  ;I \J  iCH 

[1]  CH+lpCH 

[2]  D+lipE 

[3]  i?4-“l  +  pE 

[4]  E«-((i?pl).0)\£ 

[5]  Mo,(i?+i)pi)\i: 

[6]  £[  ;  l]«»-~£,[  ;  2] 

[7]  £[  ;i?  +  2]«-~£[  ;/?  +  l] 

[8]  E-(l  0  +p£)p(  £[/>;]) 

[9]  £>(  1  0  +p£,)p(~.£’[l;]),(,£) 

[10]  1+2 

[11]  LHiJ+2 

[12]  L3 i+L2x\ElJi J]=0 

[13]  -*»Llxi(£[c/-l;  I]*£[«7;I]  )v(£[«7  +  l  ;I]*£[.7;  J]  ) 

[14]  •+Llxi(£’[*7';J  +  l];e£,[c7';J])v(£[t7';J-l]*£[t7;.r]) 

[15]  MAPZJ-liI-ll+'  ' 

[16]  -►L  2 

[17]  Ll:M4P[«7-l;I-l]«-<7ff 

[18]  L2:-»-Z,3xi(«W+l  )s(l+p£)-l 

[19]  *^L4xx  ( J4-J  +  1  );S(“l  +  pP)-l 

[20]  E+  “l  "l  +  P 

[21]  EM  1  IE 

V 


-167- 


7  CH  CONLINE  E  \  D  \  R  \  E  \  I  \J  i  CH 
Cl]  CH+lpCH 

[2]  D+lipE 

[3]  i?«-~l+p£ 

[4]  £>((7?pl),0)\F 

C  53  £<-(o,(/?+i)pi)\£ 

C  63  F[;l]-2?[;2] 

C7]  £[  ;fl+2>£[  ;/?  +  l3 

[8]  E+(l  0  +p£)p(  ,£),(  ,^CD;3) 

[9]  E+(  1  0  +p£)p(,m;3),(,£) 

[10]  1+2 

C  11]  L4:«7«-2 

[12]  L3i-+L2x\ELJ;  J]  =  0 

[13]  -*»Llxi<£[«7-i  ;J]*£[«7;I]  )v(£[«7  +  i  ;!]*£[«/;  I]  ) 

[14]  +Llx\(ElJiI+ll*ELJin)v(ElJiI-ll*ElJiIl) 

[15]  MAPiJ-1 ;  J-l  ]■«-'  1 

[16]  -*-L2 

[17]  LI iMAPZJ-1 il-ll+CH 

[18]  L2 : +L3x \(J+J+l)£(lipE)-l 

[19]  -*,L4x  x  ( I+I +l)£(*”l  +  pE,)-l 

[20]  E+  “1  "1  +£ 

[21]  E+  1  1  \E 
V 
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V  CH  CONBLANK  E\D\RiE\I\J\CH 
Cl]  CH+lpCH 
[2]  D+l*pE 
C  33  R+~lipE 
C  43  E+{(Rpl),0)\E 

C  53  0 .  (i?  +  l  )p  1 ) \£7 

[63  £C  ;l]<-£[  J2] 

C  73  EL  ;i?+2]«-£[  ;/?+l] 

[8]  E+L 1  0  +pE)p( ,E),( .ELDil) 

[9]  E+( 1  0  +pE)p(,EClil),(tE) 

[10]  1+2 

[11]  L^xJ+2 

[12]  L3  x+L2*\EL  J\ J]  =  0 

[13]  +L1*\(ELJ-1  ;!]*£[ J;I]  )v(F[tf+l;  I]*£[J;I]) 
ci4]  +Lix\(EiJii+i]*ELJin  )v(ELJii-n*EUin) 

[15]  MAPL J-l il-ll+CH 

[16]  +L2 

[17]  LliMAPZJ-liI-ll*'  ' 

[18]  L2  :  -*-L3  x  i(J-<-tl+l)^(ltpF)-l 

[19]  +L  4  x  \  (  J  +1  +  l)£(~l+p£’)-l 

[20]  E+  “l  "l  +E 

[21]  E+  1  l  +E 
7 


-169- 


V  RES+CIRCLE  XYR \R\C  ;XC \YC  iRAD  ;QiX  ;Y  -tRES 
Cl]  R+-\\pMAP 

[2]  C+-~l\pUAP 

[3]  XC+XYRll'] 

C  4]  YC+XYRl 2] 

[5]  RAD+XYRl  3] 

C  6 ]  Q+ 11*7 

C7]  X+((\C)-XC) 

[8]  Y+{ ( \R)-YC)*Q 

[9]  RES*-((Y*Y)o  .+(XxX))£RADxRAD 


V 
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V  RES+STAR  XYR\R\C\XC \YC\RADiQiE\X\Y\XYR\RES 

[1]  R+l*pMAP 

[2]  C+~l+pMAP 
C  33  XC+XYR [1] 

[4]  YC+XYR C2] 

[5]  RAD+XYR131 

[63  Q+ 11+7 

[7]  2T+2  +  3 

[8]  XH\C)-XC 

[9]  Y+((\R)~YC) 

[10]  RES+(  (Y*E)o.  +(X*E)  )<,RAD*E 

V 
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Cl] 
C  2D 
C  33 

[4] 

[5] 
C  63 

[7] 

[8] 
C  93 


V  RES+X  RECT  X ;  XR ;  XR  ;R ;  C ;  X ;  Y ;  Xl ;  X2  -tYl ;  X2  iRES 
R+lipMAP 
C+~ ltp MAP 
Xl+Xll'} 

X2+XI 2] 

yi+ni] 

X 2-*-Xl  2] 

XR+\C 


XR+\R 

RES+( (XR£X2)a(XRzX1) )°,a((XR£X2)*(XR2X1) ) 


V 
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V  BMPRINT iQl iQiIVl iIV 2 ; IVO \A ; B ; D ; P; E 

[1]  +Ll*\*/( tMAP)e(l  0) 

[2]  'USE  VMPRINT ' 

[33  -*-0 

L  Pi  :£<-l+p//i4P 

[5]  <?1«-0,1 


[6] 

[7] 

[8] 

[9] 

[10] 
[11] 
[12] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[20] 
[21] 
[  22] 

[23] 

[24] 

[25] 

[26] 


P«-  ltp MAP 
IVO+Q lp  •  |  » 

IVl+QlpQ\pQ\ (90pl0=tl0)\* 123456789* 

I72«-Glp3+(3p  *  123456  78  90') 

P«-P+( 9 Op 10  =  i 10 )\ *12  345  678  9  * 

' ENTER  TITLE • 

>4  ■*-[][] 

'ENTER  2  CHARACTERS  REPRESENTING  LEVELS  OF  GRAY' 

' POSITION  NEW  PAGE  ,  PRESS  RETURN ' 

JM3 

-►P2xi(  (p/D^P) 

*  '.A 

•*-L  3 

L  2  s  *  '  #(r(J?-p>i)*2)p  '  *  Ji4 

P3  :  '  * 

1 

'  ' ,Pp'1234567890* 

'  ' , ( P+2 )p '  * 

171,172,'  *  ,J70,fl[l+/MP],J70 
'  * , ( P+2 )p ' ” * 


7 
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7  BMPRINTliQl  1Q1IVI ; IV2; IVO'tB\C\E 
Cl]  -*»Llxi  a/(  ,MAP)e(l  0) 

[2]  'USE  VM PRINT' 

C  3]  -M) 

[4]  LI iQ+l+pMAP 

[5]  Gl-KM 

C  63  E+~l*pMAP 
C  73  IVO+Qlp  '  I  ! 

[8]  IV1 +QlpQlpQ t(90p 10=1 10 )\» 123456789* 

C  93  IV2+-Q\pQ\(  Qp  1 1234567890*  ) 

[10]  C+E+ ( 90pl0=il0)\*123456789' 

Cll]  B+',*' 

[12]  '  ' 

[13]  '  '  C 

[14]  »  » \Ep '1234567890' 

[15]  '  •#(F+2)p'„, 

[16]  IV1 % IV2 , '  * %IVOtBll+MAPl %IVO 

[17]  '  »,(£+2)p*"» 

7 
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V  CHPRINT-.Q  1  iQilVl  -.IV2\IV0-.A-,B\D-,C-,E\N\IN 
Cl]  Q+lipMAP 

C  2]  Ql+Q.l 

[3]  E+~l*pMAP 

[4]  IVO+Qlp'\  ' 

[5]  J71-«-31pClp(?‘t,(90plO::\10)\*123456789* 

[6]  IV2+QlpQHQp  *1234567890*  ) 

[7]  0£+(90pl0=il0)\'  123456789  * 

[ 8 ]  *  ENTER  TITLE ' 

[9]  X+E 

[10]  LI -.'ENTER  NUMBER  OF  DESIRED  OUTPUT  OPTION' 

[11]  ' POSITION  NEW  PAGE  AND  PRES  RETURN .OPTIONS  ARE 

[12]  *  1.  BARE  MAP ' 

[13]  '  2.  FRAMED  MAP ' 

[14]  *  3.  COORD .  FRAMED  MAP ' 

[15]  *  4.  COORD.  BARE  MAP' 

[16]  -L7xi~(//+lp[])e(  1  2  3  4) 

[17]  Itf«-(L((£-pd)*2)  )x(£>p4) 

[18]  -*-(tf-l)4>Ll,L2,L3,L4 

[19]  Ll-.INp  '  '  ;A 

[20]  *  » 

[21]  MAP 

[22]  +0 

[23]  L  2:  »  'AIN  p  *  »  ) \A 

[24]  *  * 

[25]  (£+2)p*  * 

[26]  IVO.MAP.IVO 

[27]  (£+2)p*“* 

[28]  -*-0 

[29]  L3:  (4p  9  *  )  .(Itfp  *  *  )  \A 

[30]  *  * 

[31]  '  ' »C 

[32]  *  » ,Ep *1234567890* 

[33]  *  *  . (£+2 )p  9  * 

[34]  IV1.IV2.'  ' .IVO.MAP.IVO 

[35]  '  *  , ( E+2 )p  *  *"  * 

[36]  -*-0 

[37]  L4:(3p*  *  ),(27/p*  9  );4 

[38]  *  * 

[  39 ]  *  *  .  C 

[40]  *  * .Ep *1234567890* 

[41]  *  * 

[42]  IV1.IV2.'  '.MAP 

V 
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V  VMPRINT-,NGL-.A°.I-.J-,OUTiCH',IN',T-,E'.N 

[1]  'MIN*'  ;  {l/. MAP) 

[2]  'MAX*  '  -.([/.MAP) 

[3]  19: 'ENTER  NUMBER  OF  DESIRED  OUTPUT .OPTION S  ARE  :' 

[4]  •  1.  BARE  MAP' 

[5]  '  2.  FRAMED  MAP' 

[6]  1  3.  COORD .  FRAMED  MAP ' 

[73  8  4.  COORD .  BARE  MAP' 

[8]  -*L9xi~(tf«.ip[])€(  i  2  3  4) 

[9]  +{N-1  )4>L1  ,L12,L13#L14 

[10]  L12:VMPRINT2 

[11]  -*-0 

[12]  LI 3: VMPRINT3 

[13]  *>0 

[14]  LI  4 : VMPRINT^ 

[15]  *>-0 

[16]  LI : 'ENTER  NUMBER  OF  GRAY  LEVELS' 

[17]  NGL+  loD 

[18]  'ENTER  CHARACTERS  REPRESENTING  EACH  GRAY  LEVEL  (*  BLANK)' 

[19]  CH+NGLpV\ 

[20]  'ENTER  IN  ASCENDING  ORDER  LOWER  BOUND  AND  UPPER  B.  FOR  EACH  LEVEL' 

[21]  A+(NGL. 2)p0 

[22]  1*1 

[23]  L3 :  *  LEVEL  ';I;'  REPRESENTED  BY  '  iCHUl 

[24]  4[I;]«-2pD 

[25]  +L8xiX*i 

[26]  +L1  x\/l[l  }!]<([  /  .MAP) 

[27]  -*L6xi4[l;2]^[l;l] 

[28]  *L1 

[29]  L8:-»L7xl/l[J;l3^yl[J-l  ;2] 

[30]  +L 6x\4[J;2]^4[J  » 1  ] 

[31]  L7  :->L3 ,  pD**-’  ERROR  ;  OVERLAPPING  OR  IMPROPER  BOUNDS  ;  SPECIFY  AGAIN ' 

[32]  L6  :+L3*\  (.1+1+1)  £NGL 

[33]  +L  7*\( / .A)> ([ /  . MA  P )  )v(  (L/  tA  )<([/,  MAP)  ) 

[34]  OUT+ClipMAP)pO 

[35]  'ENTER  TITLE .POSITION  NEW  PAGE. PRESS  RETURN' 

[36]  M3 

[37]  E+~ lip  MAP 

[38]  JML((^-p2,)*2)  )x(E>pT) 

[39]  INp  '  '  ;T 

[40]  1  * 

[41]  1+1 

[42]  L5: Jf-1 

[43]  Lk:OUT+OUT+(  (MAPlI  ilTzAZJ -.11)  a(MAPLI  •.  UALJ  i2l)  W 

[44]  -»-L4xi  (  J+J+l  )<.NGL 

[45]  CH+'  ' .CH 

[46]  CHll+OUTl 

[47]  0UT+( ~1 ipMAP)p 0 

[48]  CH+(1*CH) 

[49]  +LS*\U+I+l)£(lipMAP) 

7 


-176- 


7  VMPRINT2-.NGL'.A\I-.J-,OUT',CR'.IN-,T',E'.N',Q-.Ql'.E\IVO'.IV\\IV2'.C 

[I]  Q+l\pMAP 

C  23  Q1+-Q.1 

[3]  E+-~  1+pMAP 

[4]  IV0+Q1 p'|' 

[5]  I VI  -*-QlpQlpQ  +  (90pl0=\10)\* 123456789  * 

[6]  IV2+Q\pQHQp  '  1234567890*  ) 

[7]  C+E\ (9  Op 10  =  il0)\' 123456789* 

[8]  'ENTER  NUMBER  OF  GRAY  LEVELS' 

[9]  NGL+lpO 

[10]  • ENTER  CHARACTERS  REPRESENTING  EACH  GRAY  LEVEL  (*  BLANK)' 

[II]  CH+NGLpR 

[12]  'ENTER  IN  ASCENDING  ORDER  LOWER  BOUND  AND  UPPER  B .  FOR  EACH  LEVEL ' 

[13]  A+(NGL. 2)p0 

[14]  !«- 1 

[15]  £3:*  LEVEL  '  REPRESENTED  BY  '  iCHUl 

[16]  4[I;>  2pD 

[17]  ->£8xil*l 

[18]  ->£  7x\/l[l;l]<(L  /  t  P ) 

[19]  -»-L6xij4[1;2]£/1[1;1] 

[20]  -*£7 

[21]  £8  :->£7xu4[I;l]£i4[J-l  ;2] 

[22]  ->£  6x\i4[J;2]^i4[J  ;1] 

[23]  £7:->£3,pO*£i?£0£  \  OVERLAPPING  OR  IMPROPER  BOUNDS  ;  SPECIFY  AGAIN ' 

[24]  £6  :->£3xi  (  I-kT  +  1  )$NGL 

[25]  ->£7xi  (  (\ / .A)>([ / .MAP)  )  v(  (  L  /  tA  )  <  (  L  /  tMAP)  ) 

[26]  OUT+(\  lipMAP)pO 

[27]  'ENTER  TITLE .POSITION  NEW  PAGE .PRESS  RETURN' 

[28]  Ml 

[29]  IN*-(  L((£-piF)<-2)  )x(£>p2*) 

[30]  '  '.(IN p'  *); T 

[31]  •  * 

[32]  ( £+2 ) p  *  _  * 

[33]  J>1 

[34]  £5 iJ+1 

[35]  L*:OUT+OUT+( (MAPlI il^AlJ ill )  *(MAPlI  il*AiJ 2]) W 

[36]  *>£4x ;  (J-*-J  +  \)<,NGL 

[37]  CH *•  *  '  CH 

[38]  (  .IVO)\.IhCHll+OUTl  ;(  .IV0)H1 

[39]  OUT+(~lipMAP)pO 

[40]  CH+(1*CH) 

[41]  ->£5x1  (J«-J+l)£(  lip  MAP) 

[42]  (£+2 )p *  ”  * 

V 
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V  VMPRINT3 iNGLiAiliJlOUTiCHiINiTiEiNiQiQi-EilVOiIV 1 ;JT2;C 
Cl]  Q+lipMAP 

[2]  Ql+Q, 1 
C  33  E+~l*pMAP 

[4]  IVO+Qlp ' | ' 

[5]  IVl+-QlpQlpQ\  (90pl0=il0)\* 123456789  * 

C6]  IV2 +Q lpQ\(Qp  *1234567890' ) 

[7]  OF+(90pl0=il0)\,123456789  * 

[8]  *  ENTER  NUMBER  OF  GRAY  LEVELS' 

[9]  NGL+ lpD 

CIO]  * ENTER  CHARACTERS  REPRESENTING  EACH  GRAY  LEVEL  (*  BLANK)' 

C 11]  CH+NGLp 0 

C 1 2 ]  'ENTER  IN  ASCENDING  ORDER  LOWER  BOUND  AND  UPPER  B.  FOR  EACH  LEVEL * 
C 13]  A+(NGL. 2)p0 
C 1 4 ]  1+1 

C 15]  L3:' LEVEL  '  ; I;'  REPRESENTED  BY  ' iCHZH 
C  1 6 ]  i4Cl;  ]«-2pD 
C 1 7]  +LQ*\I*1 

C  1 8  ]  -+L1  *\A\.l  J 1  ]  < (  L /  .MAP ) 

C  19  ]  ■*'L6xt>lCli2]^Cl5l] 

C  20  ]  -*-L7 

C  21  ]  Z-8  :-*-L7xii4CJ»l]^i4CJ-l»2] 

C  22]  -*»L6xiXCJ;2]^/lCr;i] 

C  23  ]  L7  i-+L3  ,pD+-'  ERROR  ;  OVPPL/lPPrPG  OP  IMPROPER  BOUNDS  ;  SPECIFY  AGAIN' 

C  24 ]  L6i+L3x\(I+I+l)zNGL 

C  25  ]  -*-L7x  \(  ([  /  .A)> ([ /  .MAP  ))v((L  /  , /l  )  <  (  L  /  .MAP )  ) 

C  26]  OUT+-(  ~l+pM4P)pO 

C  27]  'ENTER  TITLE .POSITION  NEW  PAGE .PRESS  RETURN' 

C  28  ]  2VD 

C  29 ]  JP«-( L( (E-pT)*2) )*(E>pT) 

C  30  ]  (  4p  *  ').(IN  p'  ')-,T 

C  31 ]  '  » 

C  32]  *  *, C 

C  33 ]  *  • .Ep  *1234567890' 

C  34]  *  * ,(P*2)p 

C  35  ]  I«-l 
C  36]  L5:J+1 

C  37 ]  L4 : OUT+OUT+( (MAPI I ; ] ZAlJ ; 1 ] ) *(MAPl I ; ] $ALJ i 2] ) ) xj 
C  38]  •►L4xt  (  JW  +  l  )ZNGL 
C  39 ]  CH+'  ' ,CH 

C  40  ]  (,IKl)CJ],(,JV2)Ci-],  '  '  A.IVO)ZlllCHll+OUTli(  ,IVO)ZH 

C  41 ]  OUT+( ~ lip  MAP) p  0 
C  42]  CH+UiCH) 

C  43  ]  *►£  5  x  \  (  J  «-J +l)s(l+p  MA  P ) 

C  44 ]  *  '  .  (  P  +  2  )p  1 ""  * 

V 


-178- 


V  VMPRINTh  -.NGL  ',A',I  -.J  -.OUT  \CH  -.IN  \T'.E \N  '.Q-.Q1  -.E-.IVO-.  IV1-.IV2-.C 
Cl]  Q+lipMAP 

[2]  Ql^.l 

[3]  P«-”l  ipMAP 

[4]  IVO+Qlp  '  |  * 

C  5]  IVl+QlpQlpQi ( 9 Opl 0=i 10 )\» 12  3456789  * 

[6]  IV2+QlpQHQp  *1234567890*  ) 

[7]  C+Ei (9  Op 10=il0)\* 123456789* 

[8]  *  ENTER  NUMBER  OF  GRAY  LEVELS  ' 

[9]  NGL+lpU 

[10]  ' ENTER  CHARACTERS  REPRESENTING  EACH  GRAY  LEVEL  (*  BLANK)' 

[11]  CH+NGLpW 

[12]  'ENTER  IN  ASCENDING  ORDER  LOWER  BOUND  AND  UPPER  B.  FOR  EACH  LEVEL ’ 

[13]  AHNGL.  2)p0 

[14]  1+ 1 

[15]  L3:' LEVEL  '-.I-,'  REPRESENTED  BY  '  iCHUl 

[16]  4[I;]+2pD 

[17]  +LB*\I*1 

[18]  '►L  7x;>l[l;l]<(L  /  .MAP) 

[19]  -*“L6xiyl[l;2]^[l;l] 

[20]  -+L7 

[21]  L8i+L7x\AZhllsAlI-li2l 

[22]  -+L6*  ii4[J;2]^i4[J;l] 

[23]  L  7  :  -+L  3  ,  p  □«-  *  ERROR  ;  OVERLAPPING  OR  IMPROPER  BOUNDS  ;  SPECIFY  AGAIN ' 

[24]  L6:-*-L3*\  ( I+-I+1  )$NGL 

[25]  ->L7xi((r/,>l)>(r/,^P))v((L/>i4)<(L/,A//lP)) 

[26]  OUT<rCl+pMAP)p  0 

[27]  'ENTER  TITLE .POSITION  NEW  PAGE .PRESS  RETURN' 

[28]  T-HS 

[29]  IN+(  L((P“P2,)+2))x(p>p3’) 

[30]  ( 3p  *  ').(INp'  *)‘.T 

[31]  '  • 

[32]  •  '  .C 

[33]  •  • ,Pp,1234567890' 

[34]  •  1 

[35]  I«-l 

[36]  L5-.J+1 

[37]  LH:OUT+OUT+((MAPLl  i  IZALJ  ill )  a( MP[I ;  ] £d[<7 ;  2 ]  )  W 

[38]  -**£  4  x  \  ( *7  -<-*7  + 1  )  &NGL 

[39]  CU+'  ' .CH 

[40]  (.IVl)in.(.IV2)in.  '  ' -.CHll+OUT] 

[41]  OUT-*-{  -  l+  pAMP)pO 

[42]  CH+(liCH) 

[43]  -►LSxxC  I+I+l)z(lipMAP) 

V 


7  B+MREAD  NAME\R\ IND ; N 
Cl]  'MAPS'  FT IE  1 

[2]  'INDEX'  FT IE  2 

[3]  N+FFEAD  2,1 

[4]  !«-  2 

[5]  12 : IND+FPEAD  2,1 

[6]  ->£3xi(p,AMtf£)=(p,  *  ') 

[7]  -*■£  3  x  \  (  p  N  AM  E )  *  (  p  I N  D ) 

C  8  ]  -*>L1X  \  a/ N AMERIND 

C  9 ]  L3:+L2*\( I+I+1)*N+1 

CIO]  'MAP  NOT  FOUND ' 

C 1 1 ]  -*-14 

C  1 2  ]  LI :  R+-FREAD  1,1-1 
C 1 3 ]  £4 : FUNTIE  1 
C 1 4 ]  FUNTIE  2 
V 
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V  W  HISTOG  F \MAX\K\I\\I2\I3 \A\G 
Cl]  A+'  123456789* 

[2]  «XM/M0.5+F^ 

[3]  +L2*\MAX*100 

C  4]  * ENTER  NORMALIZATION  FACTOR  TO  REDUCE  MAXIMUM  TO  £  100* 

[5]  % MAXIMUM  -  * iMAX 

C  6]  ->0 

C  73  L2:+L1*\W=1 

[8]  »  !*»;*/ 

C  93  *  * 

CIO]  L1:I1«->1C1  +  (U/1X+100)] 

Cll]  l2-*-(At  *0*  )  C  1  +  (  IMAX  *10)] 

C  12]  J3-<-*  0123456789  *  C  l+MAX-(  \.MAX*10)  *10] 

C  1 3 ]  (?«-(».  *  )Ar/KHFZMAX)/\f>F)p  *  * 

C  1 4  ]  CCF+l^'D* 

C 1 5 ]  11,12,13,*  *,G 

C  16]  -+L1  x  i  (MAX+MAX- 1 )  >  0 

C 17 ]  *  *  ,  (  l+pF)p  *  *  * 

C 1 8 ]  O'  ' ,(pF)+(90pl0=il0)\»123456789 ' 

C 1 9 ]  O*  ' ,(pF)p(  *1234567890*  ) 

V 
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